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CIh' This paper presents the results of the observations of the detectors participating in the Interna- 

1~ \ tional Gravitational Event Collaboration (IGEC) from 1997 to 2000 and reviews the data analysis 

methods. The analysis is designed to search for coincident excitations in multiple detectors. The 
' data set analysed in this article covers a longer period and is more complete than that given in pre- 

\ vious reports. The current analysis is more accurate for determining the false dismissal probability 

for a time coincidence search and it optimizes the search with respect to a target amplitude and 
direction of the signal. The statistics of the accidental coincidences agrees with the model used for 
drawing the results. The observations of this IGEC search are consistent with no detection of gravi- 
tational wave burst events. A new conservative upper limit has been set on the rate of gravitational 
■ wave bursts with Fourier component H > 2 ■ W~ 21 Hz^ 1 , both for searches with and without a filter 

for the Galactic Center direction. This study confirms that the false alarm rate of the observation 
can be negligible when at least three detectors are operating simultaneously. 

PACS numbers: 0480Nn, 9585Sz 



I. INTRODUCTION 

This paper presents the results of the observations 
of the International Gravitational Event Collaboration 
(IGEC) from 1997 to 2000. We have made an extensive 
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search for burst-type gravitational waves with the largest 
network of detectors ever assembled, and report here the 
details of the search and the new upper limits achieved 
for the rate of gravitational wave burst events. 

The search for gravitational waves (gw) involves de- 
tecting the presence of a signal in the noise of the de- 
tector array. A signal must compete with the intrinsic 
noise of the detectors, and also with transient excita- 
tions (of mechanical or electromagnetic origin for exam- 
ple) which usually can not be discriminated from the 
actual gw signal. Therefore, it is not viable to per- 
form burst gw searches with a single detector. With two 
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or more detectors in simultaneous observation, the im- 
pact of local transient excitations on burst gw searches 
is significantly reduced. Moreover, the false alarm rate 
can be reliably estimated. Hence, to facilitate multi- 
detector searches for burst gw, the International Gravita- 
tional Event Collaboration (IGEC) was formed in 1997 
[32^ . This collaboration currently consists of five cryo- 
genic resonant-bar gravitational wave detectors, ALLE- 
GRO 0], AURIGA 0, EXPLORER 0, NAUTILUS 
and NIOBE Q , operating as a worldwide network. The 
members of this network exchange lists of candidate gw 
events and related information under an agreed data ex- 
change protocol. 

The target signals are transients without structure in 
the frequency range investigated. Examples of such sig- 
nals are short pulses of ~ 1 ms duration, signals showing 
a few cycles of ~ 1 ms period and signals sweeping in 
frequency across ~ 1 kHz. Possible sources are therefore 
related to compact astrophysical objects, like the coales- 
cence of neutron star and black hole binaries @, Q • 

The main method used to search for burst gw has been 
to search for an excess of coincident excitations in two 
detectors H El EH . The IGEC performed the first 
thorough search on more than two detectors on data ac- 
quired in 1997 and 1998 [l^. No claims of detection were 
made, but an improved upper limit on burst gw was set. 
In 2001 all the data acquired by IGEC members between 
1997 and 2000 were exchanged. A preliminary search was 
performed on these data |13| . 

In this article, we present the results of a comprehen- 
sive search by the IGEC for burst gravitational waves 
on the full data set. In addition to the extended ob- 
servation time, this analysis makes significant progress 
over the previous searches in the following respects: (i) 
different search thresholds are systematically tried; (ii) 
the time coincidence window is determined by the de- 
sired confidence level; (iii) a directional search strategy 
is implemented; (iv) the statistics of the estimated false 
alarms is thoroughly investigated and (v) the statistical 
methods chosen to set the confidence intervals ensure a 
given coverage |33l |. i.e. the probability that the confi- 
dence interval contains the true value. In the following 
Section, we review the IGEC operation during 1997-2000. 
The methods of the multi detector analysis are described 
in Section ITTT1 Finally, in Section IfVl we discuss the re- 
sults. In this section we pay particular attention to the 
new upper limit on the rate of detected burst gw [34[ and 
on the low level of false alarms achieved by this observa- 
tory. 



II. THE EXCHANGED DATA 

In this Section we review the 1997-2000 operation of 
the IGEC. We recall the sensitivity to burst gws of the 
participating bar detectors (see III Afl . Each detector 
group searches its data independently for gravitational 
wave candidates, or events. Then, the information ex- 



changed under IGEC is described, with particular atten- 
tion to the data validation requirements fsee 111 B|) . The 
quality of the contribution of each detector to the IGEC 
observatory is discussed in terms of observation time and 
false alarms as a function of a threshold on the ampli- 
tude of target gw signals fsec 111 C| . The statistics of the 
time series of the exchanged events show autocorrelation 
at short timescales. This clustering however disappears 
when cross-correlating different detectors fsee III D() . 



A. The detectors 

All of the currently operating resonant detectors mea- 
sure the tidal strain of a mechanically isolated cylindrical 
bar caused by impinging gws. A lighter mechanical res- 
onator is strongly coupled and tuned to the fundamental 
longitudinal mode of the bar, resulting in a system with 
two normal modes of vibration. 

A list of the main characteristics of the detectors is 
shown in Table [I] All the bars are cooled to cryogenic 
temperatures. AURIGA and NAUTILUS operate at a 
few hundred mK to further reduce contributions from 
the thermal noise. For the same reason, Niobium was 
chosen as the material for NIOBE because it has a higher 
mechanical quality factor at 4K. 

The IGEC search is focused on burst gws, which can 
be modeled as a pure Dirac (5-function excitation. The 
strength of a burst can be quantified by its Fourier am- 
plitude H 0l or amplitude, which is related to the energy 
E s deposited in the bar by 



Hn 




(1) 



where L is the bar length, M its mass, vq the mean res- 
onant frequency of the detector. 

Actually, the class of detectable signals is much wider 
than (5-functions. This search is also effective for all short 
duration signals which have an almost constant value for 
their Fourier amplitude at the detector frequencies, i.e. ~ 
700 Hz for NIOBE and - 900 Hz for the other detectors 
(see Table 0). In these cases, the gw amplitude H D is 
estimated without bias. 

The detectors were oriented to be nearly parallel to 
each other. This was done by orienting them to be per- 
pendicular to a common great circle that passes through 
or near the sites. This makes their antenna patterns co- 
herent and maximizes the probability of coincident signal 
detection between multiple detectors. 

For cylindrical bar detectors, the amplitude observed 
for a gw signal from a particular source in the sky follows 
a sin 2 function, where 9 is the angle between the long 
axis of the bar and the direction of the source. As a 
demonstration of the effectiveness of the chosen common 
orientation, we plot in Figurcnthe directional sensitivity 
with respect to the Galactic Center for the five detectors 
of the IGEC. The sensitivity to wave polarization is then 
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TABLE I: Summary of detector characteristics. The reported misalignment is the angle between the bar axis and a common 
direction. The observation time refers to the data exchanged for this 1997-2000 IGEC analysis. 



Detector 


ATT T7 1 /"1 Tl t~\ 

ALLEGRO 


A T TT> T/~1 A 

AURIGA 


EXPLORER 


"NT A TTT^TT T TO 

JNAU liLUb 


NIOBE 


Material 


A15056 


A15056 


A15056 


A15056 


Nb 


Mass [kg] 


2296 


2230 


2270 


2260 


1500 


Length [m] 


3.0 


2.9 


3.0 


3.0 


2.8 


Resonant frequencies [Hz] 


920 


930 


921 


924 


713 


895 


912 


905 


908 


694 


Temperature [K] 


4.2 


0.2 


2.6 


0.1 


5.0 


Longitude 


268°50'E 


11°56'54"E 


6°12'E 


12°40'21"E 


115°49'E 


Latitude 


30°27'N 


45°21T2"N 


46°27'N 


41°49'26"N 


-31°56N 


Azimuth 


-40° E 


44° E 


39°E 


44° E 


0° 


Misalignment [deg] 


9 


4 


2 


3 


29 


Observation time [d] 


852.5 


216.5 


551.0 


414.8 


192.6 



Directional Amplitude Response (for optimal polarization) 




5 10 15 20 



FIG. 1: Amplitude directional sensitivity of the detectors ver- 
sus UTC with respect to the Galactic Center for DAY 25 Dec 
2000. 



cos 2ip, where ^ is the polarization angle in the wavefront 
plane with respect to the projection of the bar axis. 



B. The IGEC data exchange protocol 

Each group implements independently a burst gw 
search by using an optimal filter 0], which takes into 
account the slow variations of the noise characteristics of 
the detector 0, 0, 0] . The filter estimates the Fourier 
amplitude of the burst gw. An adaptive threshold, ex- 
change threshold, is then applied to the filtered data and 
a list of the candidate events above threshold is compiled 
|3fij . The exchange threshold is typically set to ampli- 
tude signal-to- noise ratio (SNR) between 3 and 5. The 
candidate events are described by the peak amplitude of 
the optimal filter output, the time of its occurrence and 
the uncertainties on amplitude and time. 

The IGEC protocal requires each detector to report all 
time intervals of satisfactory operation. This is accom- 
plished by vetoing periods corresponding to times of lab- 
oratory activity which are known to affect the sensitivity 
of the detector, such as periods of cryogenic maintenance. 
In the case of EXPLORER and NAUTILUS, times were 
also vetoed when the noise variance exceeded a certain 



value. For AURIGA, times were vetoed when either the 
statistics of the noise was not Gaussian, or the Wiener 
filter was not properly matched to the noise [l7lll^ |. 

All this information is exchanged within IGEC under a 
common data format. This protocol was last updated in 
2000 [32]. The most relevant additions introduced by that 
update were: i) absence of any biases in the estimates 
of the time of arrival (ETA) and amplitude [35], ii) an 
estimate of the errors in ETA, iii) an upper bound of 
the systematic errors in amplitude, and iv) a continous 
measurement of the noise level and a continous record of 
the chosen exchange threshold (i.e. the threshold used 
to compile the event list). The noise level is described 
by the standard deviation and by the 3rd and 4th order 
moments of the noise distribution. 

The choice of the most suitable exchange threshold is 
left to each group, provided that two constraints are met. 
One lower bound on the exchange threshold comes from 
the requirement of having unbiased estimates of the am- 
plitude and ETA of the exchanged events. The other 
constraint is a specified criterion that limits the rate of 
exchanged events, and thus the rate of false alarms of the 
observatory. 

The search for time coincidences among events of dif- 
ferent detectors is performed by setting variable time 
windows computed from the ETA uncertainties to ensure 
a certain probability of false dismissal (see Section lllf Bj l. 
A rough estimate of the contribution of each individual 
detector to the final rate of accidental coincidences of the 
observatory can be found by considering the product of 
its average ETA standard deviation at and its average 
event rate A [13 ■ The IGEC recommendation is that the 
threshold be kept high enough so that 

ol x A = ^ l = x ° u < 0.1% (2) 

obs 

where an is the standard deviation of the arrival time for 
the ith event, n is the total number of events and T b s 
is the total observation time. This recommendation is 
an improvement of the one followed in the previous data 
exchange, i.e. to limit the rate to one hundred events 
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FIG. 2: Overview of the observation periods in 1997-2000 for each detector, as fractions of time in monthly bins. Three different 
ranges of sensitivities are considered: exchange threshold lower than 3-10 -21 Hz -1 (darker shade), included in 3 — 6T0~ 21 Hz' 1 
(middle shade) and above 6 ■ 10 -21 Hz' 1 (light shade). 




FIG. 3: Observation time of the IGEC detectors as a function 
of a threshold H t on gw amplitude in 1997-2000. The ordinate 
is the integrated time during which the detector exchange 
threshold has been lower than H t . 



values of <j t has been fractions of a second. 

Whenever environmental monitors have been operat- 
ing, the events have been checked against periods of am- 
bient disturbances prior to their exchange. If an event 
from the filtered data occurs in coincidence with an ex- 
citation observed by these monitors, it is vetoed and not 
considered a candidate gw event. For AURIGA, no veto 
based on environmental monitor has been implemented, 
yet the event has to pass a x 2 test to check its consistency 
with an impulsive mechanical excitation of the bar|l6l|. 
This test has been found to provide a good reduction 
of false alarms at high SNR. Despite all such efforts to 
remove local excitations in the detectors, most of the 
events above threshold cannot be ruled out as candidate 
gw events j22|. Partial information about vetoed events 
are tracked in the exchanged event lists for diagnostic 
purposes. [13 



per day in order to cope with the widening of the 
effective bandwidth of the detectors. 

The ETA standard deviation has been estimated by 
means of a Monte Carlo simulation for the AURIGA, EX- 
PLORER and NAUTILUS detectors HJ or by mea- 
suring the response of the detector to repeated impulsive 
excitations for the ALLEGRO detector The uncer- 
tainty on the ETA depends both on the noise level and on 
the timing accuracy of the filtered data. Therefore, the 
behavior of at is significantly different for the different 
detectors, though it decreases in all of them as SNR in- 
creases. It was also found to vary with time in AURIGA, 
EXPLORER and NAUTILUS, following the variations of 
the effective bandwidth of the detector [20l l2lj . Typical 



C. Data overview 

The duty cycle of the network has been previously re- 
ported in detail 0, 0] . A graphical representation of 
the on-off times of the individual detectors between 1997 
and 2000 is shown in FIG. |21 During this 4 year pe- 
riod, there were 1319 days when at least 1 detector was 
operating, 707 days with at least 2 detectors in simulta- 
neous operation, 173 days with at least 3 detectors and 
26 days with at least 4 detectors. This time coverage 
is a consequence of the non optimal overlap among the 
operating times of the single detectors and of their duty 
cycles, which were included within the 0.9 achieved by 
ALLEGRO and the 0.3 of AURIGA. 
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FIG. 4: Cumulative distribution function of the amplitude of the exchanged events: the ordinate is the fraction of events whose 
amplitude is greater than the value in abscissa. The amplitude is normalized to the corresponding exchange threshold of the 
detector. 
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FIG. 5: Average rate of exchanged events per each UTC hour 
of the day. From top to bottom: ALLEGRO, AURIGA, EX- 
PLORER, NAUTILUS and NIOBE. 



The cumulative observation time of each detector ver- 
sus gw amplitude is shown in FIG. The exchange 
thresholds of the detectors fluctuated significantly in time 
for all detectors but ALLEGRO, following the non sta- 
tionary behavior of the noise. 

The typical range of the exchanged thresholds was 
2-^6- 1CP 21 Hz" 1 . Neglecting the directional sensitivity 
of the detector, this range corresponds to a 1ms burst 
generated by 0.02 4- 0.2M Q solar masses converted in gw 
with isotropic emission at a distance of lOfcpc. With re- 
spect to the first IGEC data exchange [HI, EXPLORER 
and NAUTILUS take the opportunity of the relaxed rec- 
ommendation on the average rate of exchanged events 
(Eq. [2J) , and lower the exchange threshold (from SNR ss 
5 to about 4.5 ). The data of the other detectors keep the 
previous exchange thresholds (SNR « 3 for ALLEGRO 
and NIOBE, SNR « 5 for AURIGA). The mean rate 
of events exchanged by EXPLORER and NAUTILUS 
is about 5 times greater than in the previous exchange, 
while for the other detectors it remains at the same level. 

The amplitude distributions of the exchanged events 
are shown in FIG. At least two distinct regimes are 
observed in all detectors, though showing a large vari- 
ability among them: a steep roll-off close to the thresh- 



old, and an additional tail dominating at SNR greater 
than ~ 10. The exchanged events are mostly generated 
by non-Gaussian noise sources, which are not currently 
modeled. 

The average number of events observed in each hour 
(Universal Time) is shown for the five detectors in FIG.|3] 
All detectors show an increase in the number of events 
observed at certain hours of the day. For AURIGA, EX- 
PLORER and NAUTILUS, higher event rates are ob- 
served between 5 and 18 hours local time, while for 
NIOBE between 7 and 16 hours local time. Such behavior 
is consistent with a correlation with human activity. The 
event rate for ALLEGRO is almost constant and shows a 
small increase between 18 and 23 hours local time. Note 
that since AURIGA, EXPLORER and NAUTILUS are 
in the same time zone, the rise and fall of event rates is 
almost synchronized. 

The average rates of the events exchanged by each de- 
tector are shown in FIG. HO as a function of an absolute 
search threshold on the gw amplitude. The search thresh- 
old is applied to each event list before performing the co- 
incidence search described in Section Till Bl For a given 
search threshold value H t , the rate is computed dividing 
the number of events exceeding H t by the observation 
time during which the exchange threshold of the detector 
has been lower than H t . This procedure is consistent with 
the data selection we applied in the multi detector anal- 
ysis, as described in Section fill Al In general, the event 
rates decrease as the search threshold increases because 
the number of selected events decreases and the selected 
observation time increases. However, this can happen in 
a non monotonic way due to the non-stationary behavior 
of the noise performances of the detectors. In fact, as the 
search threshold increases, the selected observation time 
can extend to additional periods of operation character- 
ized by worse sensitivity. These periods contribute with 
a much higher instantaneous rate of events, given that 
most events appear very close to the exchange thresh- 
old (see FIG. 0}. Therefore, the mean event rate may 
increase at some higher search threshold. 

The mean timing uncertainty, o7, of the exchanged 
events is shown in FIG. for each detector as a function 
of the search threshold. The value of 07 is dominated 
by the selected events which are closest to the exchange 
threshold of the detector. As the threshold increases, 
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FIG. 6: Average rate of the events exchanged by the detectors 
as a function of a threshold Ht on gw amplitude, as described 
in the text. The minimum threshold plotted for each detector 
corresponds to 10 days of observation time. 



the mean timing uncertainty decreases down to the limit 
given by the timing calibration or resolution of the detec- 
tors, which is in the range 1 -=-80 ms. Similarly to FIG. El 
at is a non-monotonic function of the search threshold 
due to noise being non-stationary. For instance, the peak 
shown by AURIGA data around 10~ 20 -ffz -1 is due to a 
one week period of operation of the detector with re- 
duced sensitivity and reduced effective bandwidth. For 
the NIOBE data a conservative estimate of at was pro- 
vided, independent from the event amplitude. 

FIG. and FIG. allow the reader to compare the 
contributions of each detector to the false alarm rate of 
the observatory as a function of the selected threshold H t 
on gw amplitude. In this respect, the quality of the de- 
tector is given by the product of the mean event rate and 
the mean uncertainty of the ETA, as already discussed in 
connection with eq. |2 The cleanest detectors have been 
ALLEGRO and AURIGA in the investigated range of gw 
amplitudes. To be precise, the actual false alarm rate of 
the observatory is not directly related to the mean event 
rates. Rather it is time by time proportional to the prod- 
uct of the instantaneous event rates of the participating 
detectors |2l|. Moreover, the dependence of the actual 
false alarm rate on the uncertainties of the ETA is also 
not simple: since the larger timing uncertainties dom- 
inate the time window used for the coincidence search 
(see Section Till B|l . the role of the detectors showing the 
worst timing performance is enhanced. 



D. Event time series statistics 

In a time coincidence search, the statistics of the esti- 
mated time of arrival of the events plays a fundamental 
role. In case the event times are random, a Poisson point 
process would fit the data. Actually, the data do not re- 
produce a homogeneous (i.e. stationary) point process, 
due to the changing performance of the detector and to 



FIG. 7: Average standard deviation of the ETA of exchanged 
events, at, as a function of a threshold Ht on gw amplitude, 
as described in the text. The minimum threshold plotted for 
each detector corresponds to 10 days of observation time. For 
NIOBE only a Is upper limit for the time standard deviation 
is available. 

the statistics of outliers. For our purposes, it is enough to 
check timescales below ~ 1 hour, because they are crit- 
ical for the method implemented to estimate the noise 
background (see Section lill D|) . 

In order to investigate the ETAs statistics of IGEC ex- 
changed data we calculate their correlation histograms 
(or correlograms) as shown in FIG. Correlograms are 
histograms of the time lags between ETAs. For a Poisson 
process one expects the histogram to be flat -i.e. without 
preferred time delays between events (see appendix lAlfor 
details). It turns out that the ETAs show relevant auto- 
correlation at small timescales, down to a few seconds. 
A clustering of the event time series for all the five IGEC 
detectors is evident. The clustering disappears as soon 
as one looks to the cross-correlation properties. This is 
very relevant for the next phase (i.e. coincidence search), 
because if this were not the case, the output of the em- 
pirical method we use to estimate the background would 
be biased (see Section Hit D|) . 



III. MULTIPLE DETECTOR ANALYSIS 

In this Section we present the methods implemented 
by IGEC to analyze the exchanged data. The analysis is 
based on a time coincidence search among gw candidates, 
or events, of different detectors (see 1111 B")l . The time co- 
incidence window is varied to get the desired maximum 
probability of false dismissal. There is also a test of the 
compatibility of the signal amplitude, as estimated by 
the different detectors (see IIII CJI . Prior to the coinci- 
dence search, we apply a data selection procedure which 
limits the search to burst gws exceeding a specified search 
threshold fsee llll A"j) . A directional search strategy is im- 
plemented as well. The IGEC analysis is then performed 
as a function of the search threshold value and of the gw 
direction. The main advantages of this procedure are a 
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FIG. 8: Correlogram matrix for the estimated time of arrival of all exchanged events in the years 1997-2000; all plots are 
histograms of the time lags between events in units of counts/sec (ordinate) vs seconds (abscissa). All five auto-correlation 
histograms of the single detectors (diagonal plots) show correlation at short timescales, and in three cases these structures 
dominate by far over the uniform distribution expected in case of a Poisson process. Despite this, the cross-correlograms 
(off-diagonal plots) show no sign of residual correlation. The NA-NI cross-correlogram is empty because these two detectors 
had no overlap in their operative time. 



reduction in the false alarm rate and the control of the 
maximum probability of false dismissal of any burst gw 
exceeding the search threshold. The accidental coinci- 
dence background has been estimated by applying the 
same analysis procedures on many data sets obtained by 
shifting the time of the real data (see IIII Dfl . The first 
relevant result of this analysis is that the number of esti- 
mated accidental coincidences turns out to be a Poisson 
random variable. Finally, in IIII El we describe the statis- 
tical method used to set the confidence intervals on the 
number of gw bursts detected by the observatory. This 
method is unified and ensures the desired coverage of the 
resulting confidence intervals. 



A. Data selection 

Before searching time coincidences, we apply a data 
selection procedure which limits the search to burst gws 



exceeding a specified amplitude. 

The first step is to specify an absolute threshold H t , 
or search threshold, for the gw amplitude estimates. This 
threshold is common to all detectors and sets the lower 
bound on the amplitude of the target gw population. The 
following multiple detector analysis is then repeated sys- 
tematically for different threshold values. 

The second step is to exclude from the observation 
time of each detector all time periods where the exchange 
threshold is above the chosen search threshold. The mo- 
tivation is to limit the false dismissal probability of any 
burst gw of amplitude greater than the selected thresh- 
old. With this selection, the detection efficiency for burst 
gws exceeding H t is at least 0.5 in any detector. 

The third step is to exclude candidate events that are 
below the search threshold. This gives a significant re- 
duction of the false alarms, while preserving the same 
minimum detection efficiency of the previous step. In 
fact, the exclusion of the lower amplitude events cuts 
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down both the rates and the time uncertainties of the 
events of each detector (see FIG. El and FIG. [7}. 

This procedure differs from what have been previously 
done in the field j^]. The data selection is illustrated 
in FIG. El with a sample of AURIGA data. In general, a 
higher search threshold H t allows new portions of the ob- 
servation time to be considered, thus increasing the effec- 
tive observation time. Moreover, the rate of the selected 
events will strongly depend on the value of the detector 
exchange threshold with respect to H t . In particular, 
higher rates of events are favored whenever the exchange 
threshold approaches and crosses the search threshold. 

The described data selection is suitable for a blind 
search over the sky, irrespective of the source location. 
In fact, a burst gw is seen at each detector with the same 
amplitude, given that the antenna patterns are almost 
coherent. 

It is also possible to implement a directional search 
strategy to optimize the search for a specific gw direc- 
tion. This has been accomplished by modulating the ex- 
changed data with the directional sensitivity of the detec- 
tors. Specifically, all exchanged amplitudes (event am- 
plitudes, exchange thresholds, etc.) are divided by the 
time-dependent angular attenuation factor for the spe- 
cific direction in the sky (see FIG. In this way, all 
amplitudes are given in terms of a burst gw propagating 
from the selected direction. FIG. 1 101 shows the search for 
burst gws from the Galactic Center direction on the same 
data set of FIG. El Then, the rest of the selection proce- 
dure above is applied. The effective observation time is 
reduced with respect to that obtained for a blind search 
over the sky at the same threshold H t because the pe- 
riods when the detector is not favourably aligned with 
the source are removed. Additionally, the set of selected 
events is generally different. 

As a result, the background noise is reduced while the 
detection efficiency for the selected direction is preserved. 
It is worth noticing that the angular selectivity of any 
directional search is quite poor, due to the broadness of 
the directional sensitivities of the detectors. For instance, 
when the detector is optimally oriented with respect to 
the chosen direction, any source within ±20deg from it 
is seen with at most 11% attenuation. 

The modulation by the directional sensitivity corre- 
lates the amplitudes of events and exchange thresholds 
among different detectors in a much more significant way 
than any other observed daily effect. This modulation 
produces new cuts to the observation time and related 
clusters of events which are almost synchronized in dif- 
ferent detectors. As a consequence, the probability of 
coincidences is enhanced at the edges of the time spans 
of common observation. The effects of this when esti- 
mating the rate of accidental coincidences is discussed in 
Sec. lrriDl 

We note that the false dismissal contributed by this 
data selection is at most 50% per each detector for burst 
gws of amplitude exceeding the chosen search threshold 
H t . The overall false dismissal of this multiple detector 
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FIG. 9: An example of the selection on AURIGA data for a 
search not optimized for a specific gw direction. The continu- 
ous gray line shows the amplitude of the exchange threshold, 
the dots represent the exchanged events vs. time. In the dark 
gray period no data was exchanged. The light gray shaded 
area shows the amplitude-time regions excluded by the data 
selection at a search threshold Ht = 3-10~ 21 Hz 1 ■ The obser- 
vation times with exchange threshold > 3-10 -21 Hz -1 are now 
excluded as well as the events with amplitude < 3-1CP 21 Hz -1 . 
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FIG. 10: The same data of FIG.|5]but now specialized to a 
search for gravitational wave bursts from the Galactic Center. 
The exchange threshold (gray line) and events (dots) have 
been divided by the amplitude directional sensitivity of the 
detector to the Galactic Center. The resulting amplitudes are 
in terms of a burst gw from the Galactic Center. The light 
gray shaded area shows the amplitude-time regions excluded 
by the data selection at a search threshold Ht = 3-10~ 21 Hz" 1 . 



analysis will be further increased by the next two steps, 
namely the time coincidence search and the amplitude 
consistency check (see the following Sections). However, 
the overall false dismissal for a burst gw above H t is un- 
der control at least in a conservative sense. The optimal 
range of H t values that would enhance the chances of a 
gw detection will depend also on the (unknown) ampli- 
tude distribution and rate of the burst gws. 
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B. Time coincidence search 

Two events from different detectors are defined to be 
in coincidence if their estimated times of arrival U and tj 
are compatible within their variances of. and of. : 

\U -tj\< At i:j = (k^jcrl+a 2 tj + M M ax) (3) 

where Atij is the coincidence window and AtMAX is 
the maximum expected light travel time between the de- 
tectors. At^ is computed according to the desired proba- 
bility of false dismissal by setting k through a Bicnaymc- 
Tchebyscheff inequality. 

Specifically, the probability to miss a coincidence be- 
cause of the uncertainties in the ETA is 

FD = P{\ti-tj\ > At i:j } 

which in turn is upper bounded by 

FD < P{\U -tj\> ky/al+a?.} < ^ = P T (4) 

where Pt is the maximum false dismissal probability 
chosen for the coincidence search (see Appendix [BJ). The 
resulting coincidence window changes for each couple of 
events according to their ETA variances. 

This procedure ensures that the required false dis- 
missal probability is met regardless of the distribution 
of the estimated arrival time uncertainties. In fact, it is 
not possible in general to approximate the ETA statistics 
as Gaussian, due to either the intrinsic narrow bandwidth 
of the detectors or their limited time resolution. For in- 
stance, in the AURIGA data the time accuracy has been 
1 ms and the time uncertainty distribution is multi-modal 
[23| at the SNR of interest. The ETA standard deviation 
is then dependent on the SNR of the event. We note that 
a t ~3> At max is the standard condition in these IGEC 
data. For comparison, the previous IGEC preliminary 
analysis 01 nas been performed with a fixed value for 
the coincidence time window and therefore only a vague 
indication of the false dismissal was possible. 

The implemented coincidence search allows one event 
to be in coincidence with more than one event in the other 
detector. A coincidence in more than two detectors has 
to satisfy eq.|3]for all combinations of detector pairs, and 
the resulting conservative false dismissal has to take into 
account the number of such required conditions. The 
same time coincidence search algorithm has been applied 
both to actually search for burst gws and to estimate the 
corresponding number of accidentals (see Section IIII D|) . 

The choice of the conservative false dismissal Pt due 
to the coincidence search can be optimized (see Ap- 
pendix [5]) . It turns out that Pt should be set between 
30% and 5% to achieve a satisfactory balance between 
false alarm and false dismissal probabilities in a two-fold 
coincidence search. 



C. Amplitude consistency check 

The estimated amplitude of a burst gw is affected by 
the systematic and statistical uncertainties of the detec- 
tor. Hence, once a set of events are found to be in coin- 
cidence, one can test if the differences of their estimated 
amplitudes are consistent with zero within a chosen con- 
fidence level. The goal is to lower the false alarms by 
removing the accidental coincidences whose amplitudes 
are not consistent. The test we implemented is similar 
to the condition described above to define a time coinci- 
dence Hlj . It takes into account both the variance o\. 

and the 4th central moment of the estimated ampli- 
tude of the events Ai, which are included in the current 
IGEC exchange protocol. 

Two events from different detectors have consistent 
amplitudes if 

\Ai -Aj\< AAij + AA t + AAj (5) 

where 

and AAi are the systematic amplitude calibration er- 
rors of the detectors. Pa is the required conservative false 
dismissal of the test for a coincidence which corresponds 
to a burst gw. The two alternative terms in curly brack- 
ets come from the Bienayme inequality of 2nd and 4th 
order respectively applied to the random variable Ai—Aj . 
The more stringent of them is chosen time by time (see 
Appendix [5J| • 

The general Bienayme inequality is used because the 
amplitude noise distributions of the detectors were not 
Gaussian nor modeled for a significant fraction of the 
observation time. This test is conservative and provides 
a less stringent removal of false alarms at low SNR am- 
plitudes with respect to tests based on some modeled 
statistics. Instead, at high SNR amplitudes the system- 
atic calibration errors, 10% for all detectors, dominate 
over AA^. 

The efficiency of the false alarm rejection of this am- 
plitude consistency test was found to be strongly depen- 
dent on the search threshold. The false alarms were 
significantly reduced only at high thresholds, H t > I • 
10 _20 i?z _1 . This is due to two concurrent facts. First 
of all, the implemented data selection forces most of the 
events to have similar amplitudes, since they are con- 
strained from below by the imposed search threshold and 
from above by the steep slope of the amplitude distribu- 
tion (Fig. 0}. Therefore, the data preprocessing itself 
provides an implicit rejection of most of the events in 
coincidence which show non-consistent amplitudes. Sec- 
ond, the efficiency of the test is greater at high SNR 
amplitudes, because there the amplitude differences of 
accidental events can be larger in terms of standard de- 
viations. 
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As a result, the test turned out to be convenient only 
at high search thresholds and with values of conserva- 
tive false dismissal Pa < 30%. On the contrary, at low 
thresholds the implementation of the test is disadvanta- 
geous because the increase in the false dismissal is not 
balanced by a sufficient false alarm reduction. For the 
results described in this paper, the application of this 
amplitude consistency test did not add new significant 
information. 



D. Background estimation 

In order to assess the statistical significance of the 
number of detected coincidences, a reliable estimate of 
the background -i.e. the number of coincidences that 
have been found by chance- is needed. The ideal ap- 
proach would be to obtain new independent data sets 
from a population with the same statistics but with all 
gw sources switched off, to repeat the search procedure, 
and then to compare the statistics of the found coinci- 
dences with those obtained from the original sample. 

If the ergodic hypothesis applies to our data, a good 
method to create an independent data sample is to per- 
form a relative translation of the time coordinate of the 
data exchanged by different detectors. This operation 
preserves the statistics of the single event list (average 
number of events, instantaneous rate fluctuations, auto- 
correlation of event times, ...). The coincidence counts 
found on time shifted data sets are independent as long 
as the applied time delays are longer than the maximum 
time window used in the coincidence search. It is rea- 
sonable to assume that the gw events are a negligible 
fraction of the total number of events at each detec- 
tor and therefore the gw events will not affect the co- 
incidences within the shifted data sets. This empirical 
method has been widely used in the field (see for instance, 
II HE IH 0) We remark that this framework allows 
us to only estimate the level of the coincidence back- 
ground, but not to distinguish between coincidences due 
to gw and those due to other common sources. 

There are a number of technical subtleties to address 
in order to give a complete description of the way this 
method has been practically implemented in this analy- 
sis. One issue regards the common observation time of 
each time shifted configuration, i.e. the collective length 
of the time spans after the data selection phase. In fact, 
in case of significant changes of the shifted observation 
time one should consider the rates rather than the counts 
of accidental coincidences, but then, the statistics of this 
random variable is no longer Poissonian. 

Another crucial issue is to check the statistics of the 
background, especially in IGEC data since the statistics 
of the event time series of each detector shows evidence 
for auto-correlation at small timescales (see FIG.|5J). We 
require that the statistics of the coincidence background 
estimates be Poissonian and stationary for any applied 
time shift, as is expected if the coincidence time series 




coincidences # 

FIG. 11: Sample of the histogram of the background coin- 
cidences obtained by the time shift method with a superim- 
posed Poisson fit. The plot shows the shifted coincidences be- 
tween EXPLORER and NAUTILUS searched with P T = 5%, 
Pa = and search threshold of 5.62 ■ 1CP 21 Hz' 1 without op- 
timizing for a gw direction. The one tail \ 2 probability of the 
sample is 0.7. The x 2 test is performed only on the histogram 
bins with at least 10 counts. The maximum time shift reached 
±5000s in 10s steps. 



can be modeled by a Poisson point process. In fact, this 
is the model we use in Section IIII El to estimate the statis- 
tical significance of the found coincidences with respect to 
the background. Note that the instantaneous coincidence 
rate may be also time- varying (i.e. non- homogeneous), 
and still the coincidence counts in a fixed time span of 
a shifted configuration would be a sample of a Poisson 
random variable. Independence (i.e. randomness) of suc- 
cessive coincidence times is the key to guarantee this re- 
sult. This requirement is met when the event lists are 
not cross-correlated within timescales up to the maxi- 
mum applied time shift. 

When performing a search optimized to a specific di- 
rection (Section IIII Af) . we pointed out the additional 
problem related to the appearance of clusters of events 
at times correlated in different detectors. In this case 
the key for a successful estimation of the background has 
been to time shift the data before applying the amplitude 
modulation. 

In this analysis, the number of tested time-shifted con- 
figurations for a pair of detectors is N s ~ 1000. This 
number is limited by two requirements. The first is to 
keep the shift step larger than the longer time window 
used for coincidence search. On the other hand, the max- 
imum delay has to be small enough to keep stationary the 
random variable number of accidental coincidences. The 
maximum time shift has been limited to ~ 5 • 10 3 seconds 
and the shift steps has been between 5 and 15 seconds. 

The easiest way to check the Poisson statistics and 
the independence of the background samples is to cre- 
ate an histogram of the number of coincidences found at 
each time shift and fit it to a Poisson density function. 
FIG. ^] shows the agreement of the background statis- 
tics to the model for a sample configuration. To test the 
goodness of the Poisson fits we applied chi-square tests 
whenever the histograms of the background counts were 
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FIG. 12: Histogram of the significance level of the goodness- 
of-the-fit test for the Poisson model of the background es- 
timates (one tail \ 2 probabilities). The \ 2 test has been 
performed on all the configurations of the observatory which 
ensured at least one degree of freedom (see FIG. \TT\ . The 
plotted data include configurations of different detectors with 
different search threshold values and same conservative false 
dismissal of the coincidence search (Pt = 5%, Pa = 0), with- 
out optimizing for a gw direction. The maximum time shift 
reached 5000s in 10s steps. The histogram is well described 
by a uniform distribution, so we can conclude that the x 2 
tests give results consistent with a general agreement of the 
background statistics to the Poisson model. 



sufficiently populated. The probability of getting a chi- 
square greater than the observed value has to be a sample 
of a uniform density between and 1 if the performed fit 
is good. This has been actually confirmed, as shown by 
the histogram of the resulting one-tail chi-square proba- 
bilities shown in FIG^l therefore, we can conclude that 
no deviations of the background estimates from the ex- 
pected Poisson statistics are observed. 

We remark that this is a quite relevant result, since 
in our case the single event lists showed evident auto- 
correlation at short timescales (see FIG.|3J) and since the 
gw search from a specific direction brings in a significant 
correlation of the selected event rates in different detec- 
tors, as described in Section lill Al As a consequence, we 
have been able to use all the available data to set con- 
fidence intervals on the detected gw with the procedure 
described in the following Section. 



scribed in the previous subsection. Having observed a 
number of actual coincidences N c , the likelihood func- 
tion of the average number of detected gw, AT A , is: 



*(JV A ; N c , N b ) = 
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Then the confidence interval is defined by integrat- 
ing the likelihood over the smaller domain [N- mi -j- A^ sup ] 
which ensures that the integral amounts to a specified 
value I. 



I = 



£(N)dN 



JVa, 



£(N)dN 



(8) 



The intervals defined as such have nice properties [25| . 
First, they are naturally bound to the physical domain. 
Moreover, they include the most likely estimate of the 
number of detected gw, which is zero when the expected 
background exceeds the found number of coincidences. 
Such most likely confidence intervals are also the most 
credible in the Baycsian framework assuming a uniform 
prior for AT A > 0. However, when the value for the like- 
lihood integral I is properly chosen, the resulting inter- 
val has also a well-defined minimum frequentist coverage. 
From numerical computations it turns out that / = 0.94 
and / = 0.97 guarantee the coverage to be at least 0.90 
and 0.95 respectively. By taking into account the ef- 
fective observation time, the confidence intervals can be 
expressed in terms of the Poisson rate of detected gws. 

The overcoverage ensured by this unified method is 
quite significant for true gw rates much less than the 
background rate. In particular, in case the true value 
is exactly zero (null hypothesis) the complement of the 
coverage to unity can be interpreted as the false detection 
probability, and it is lower than ~ 4.5% and ~ 2.3% when 
the conservative coverage is respectively 0.90 and 0.95. 



E. Setting confidence intervals on detected GW 

The found number of coincidences and the expected 
background are compared under the hypothesis of a su- 
perimposed homogeneous Poisson rate of detected gw sig- 
nals. The results are then expressed as confidence inter- 
vals containing the true rate of detected gw with a given 
probability, i.e. assuring a given coverage. The procedure 
we adopt to set confidence intervals follows the track of 
previously reported unified methods pil l25j . The results 
will be interpreted cither as upper limits or as evidence 
for detection in case the confidence intervals include or 
not the null result. 

We assume a Poisson model for the coincidence back- 
ground, whose expected number N& is estimated as de- 



IV. RESULTS 

The IGEC observations have been analyzed both by 
performing a blind search over the sky (without selecting 
a specific gw direction) and by optimizing the search for 
burst gws from the Galactic Center direction. In the first 
case the results refer to the amplitude component of the 
burst gw along the detector axes. In the second case the 
results are given in terms of the amplitude of a burst 
gw from the Galactic Center. In both cases, only the 
polarization component along the bar axis is considered. 

We define the operating time of a particular configu- 
ration of detectors to be the subset of the network op- 
eration periods when only the detectors of this config- 
uration are simultaneously operative. The main advan- 
tage of this procedure is that the results from different 
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configurations are then automatically independent, since 
they refer to disjoint observation times. We remark that 
the operating times of the configurations depend on the 
search threshold (see Section IITi"£|l . Within IGEC, 18 
different configurations of detectors have been operating 
during 1997-2000: 9 pairs, 7 triples and 2 four-fold con- 
figurations. The multi detector data analysis has been 
performed separately for each configuration as a function 
of the search threshold H t in the range 2 — 50-10 -21 -ff,z -1 . 

In order to synthesize the overall result of the obser- 
vatory at each threshold value, we sum the observation 
times, the coincidence counts and the backgrounds over 
all the configurations at the same H t . A confidence inter- 
val of the whole network is then re-computed accordingly 
for each H t . 

We report a synthesis of the results (number of coinci- 
dences, background, observation time, confidence inter- 
val on detected gws) for each configuration of detectors 
and investigated threshold value in Appendix [U] The 
results depend also on the choice of the parameters of 
the analysis, namely the false dismissal on the time co- 
incidence search, the false dismissal on amplitude consis- 
tency of events in coincidence, the selection of a direction 
in the sky and the required conservative probability, or 
coverage, of the confidence intervals. 

The results obtained at each threshold value are cu- 
mulative, i.e. they apply to detected burst gw whose 
amplitudes are > H t . This is a consequence of the data 
selection described in Section fill Al In particular, an up- 
per limit set at some threshold H t is valid with the same 
conservative coverage for any higher threshold. 



A. No statistical evidence for detected GW 

The overall results over the entire time span 1997-2000 
are well in agreement with the estimated background. In 
fact, the resulting confidence intervals on the number of 
detected gw signals include the null result in almost all 
the many trials performed (see Appendix Q. Only a few 
2-fold configurations give gw detections for some specific 
values of the parameters of the analysis. We will show, 
however, that the relative frequency of these cases is well 
accounted for by the probability of false alarm, i.e. of 
getting by chance a detection in case no gw were present 
in the data. 

The probability that a confidence interval may fail to 
include the null result in case no gw are present in the 
data has been numerically computed by simulating the 
procedure to set confidence intervals (see Section Till Efl . 
Due to the overcoverage of these confidence intervals, this 
probability is always smaller than the maximum false dis- 
missal associated to the interval, i.e. 1 — coverage, es- 
pecially at low background levels. Specifically, when the 
accidental coincidences are Nj, > 1, the expected num- 
ber of false alarms oscillates around ~ 3% and ~ 1.5%, 
respectively for maximum false dismissal values of 10% 
and 5%. This means that false alarms of the order of 



one every 30 (resp. 70) independent trials are the rule in 
case of N& > 1. Instead, in the limit of low background, 
i.e. N(, < 0.01, the false alarm probability turns out to 
be Nfc regardless of the required coverage. Therefore no 
false alarms are expected for the configurations showing 
low enough background levels. These predictions on false 
alarms obtained by numerical simulations have also been 
confirmed by an independent empirical method (see Ap- 
pendix Q ■ 

No statistical evidence for detected gw has been found. 
In fact, on one hand the null result is always included in 
the confidence intervals corresponding to low background 
levels, such as those related to triple and four-fold config- 
urations of detectors. On the other side, the total number 
of detections is consistent with the total number of ex- 
pected false alarms. For the many trials reported in Ap- 
pendix El the expected {found} total numbers of false 
alarms for the Galactic Center search arc 3.5 {2} and 1.7 
{0} for 0.9 and 0.95 coverage respectively. For the blind 
search over the sky these numbers are 2.5 {4} and 1.4 
{2}- 

With at least three IGEC detectors in simultaneous 
operation the false alarms are extremely rare even at 
search thresholds H t close to the exchange thresholds of 
the single detectors. So, even after many years of ob- 
servation time, such configurations would allow to easily 
identify any detected gw. Instead, in the case of two- 
fold coincidence search, the false alarms populate the 
achieved observation time up to high search thresholds, 
H t < 1Q- 20 Hz- 1 . 



B. Upper limit on the rate of detected GW 

The resulting IGEC observations can be synthesized 
by an upper limit on the rate of detected burst gw, mod- 
eled as a Poisson point process with constant rate. This 
upper limit as a function of the signal search threshold 
is given in FIG. ^5] and ^] for a search not optimized 
for any specific direction and optimized for the Galac- 
tic Center direction respectively. These upper limits are 
based on the confidence intervals of the network, com- 
puted from the sums of the data (observation times, coin- 
cidence counts and background) of all configurations per 
each threshold value. The plotted upper bounds have a 
probability of at least 95% to be greater than the actual 
gw rate value, so that the dashed region is excluded with 
the same confidence. The lower limits of the confidence 
intervals are all at null gw rate but one per each type 
of search (Galactic Center and blind). These detections 
can however be explained as expected false alarms (see 
previous Section and Appendix [UJ and do not affect the 
upper limits shown in the Figures. 

The upper bounds set by the network show in a few 
cases higher gw rates at higher thresholds than at lower 
thresholds (sec Appendix Q. This happens for the same 
reason why Fig. and Fig. may show increasing event 
rates and timing errors as H t increases. In these cases, we 
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FIG. 13: Upper limit for the rate of the burst gws detected by 
the observatory as a function of the amplitude search thresh- 
old H t . No directional search has been applied (i.e. blind 
search over the sky) . The Fourier amplitude in abscissa refers 
to the component of a burst gw along the axes of the detec- 
tors. The dashed region above the continuous line is excluded 
with at least 95% probability. The uncertainties on the esti- 
mated background affect neglibly the plotted curve. The time 
window for the coincidence search has been selected to limit 
the related false dismissal probability to at most 5% and no 
test on the amplitude consistency between events has been 
applied. The maximum amplitude systematic error related 
to the calibrations of the detectors is shown as the error bar 
parallel to the abscissa. 



chose to consider the most stringent upper bound value, 
on the basis of the fact that an upper limit computed at 
some threshold is also valid for any higher value of the 
threshold - as already remarked. This choice introduces a 
marginal bias on the stated coverage close to the plotted 
upper bounds. 

The results of the blind search, FIG. show a flat 
upper limit on the gw rate at high search thresholds, 
H t > 10 _20 i/z _1 , where no coincidences have been 
found. This rate is determined by the total observation 
time of the network at those search thresholds, T b s and 
takes into account the conservative false dismissal of the 
time coincidence search. Specifically, the rate is given 
by F /{T b s CL C ), where CL C is the confidence level of 
the time coincidence search procedure, CL C = 0.95 for 
the data plotted, and the factor F depends on the re- 
quired coverage, F = 3.6 for 0.95 probability. The confi- 
dence intervals on the rate widen at intermediate search 
thresholds mainly because of the presence of accidental 
coincidences. At Ht < 3 ■ 10~ 21 the upper limit sharply 
increases because of the corresponding cut off on the ob- 
servation time. The results referred to the Galactic Cen- 
ter direction, FIG. El show a similar behavior, smoothed 
by the effect of the modulation of the directional sensitiv- 
ity of the detectors. The relevant data are also tabulated 
in Appendix C. 



FIG. 14: Same as FIG. 1131 but specialized for the rate of 
burst gws from the Galactic Center direction. The Fourier 
amplitude in abscissa refers to the amplitude of a burst gw 
from the Galactic Center direction with optimal polarization 
with respect to the detectors. 



C. Final remarks 

The reported results refer to the gw detected by the ob- 
servatory, that is to say to the gw present as coincidences 
in the IGEC exchanged data during the 1997-2000 obser- 
vation time. In order to extend these results to the flux 
of gw crossing the earth, one should take into account 
the actual efficiency of detection, which also depends on 
the specific model of the gw source, in particular on the 
gw amplitude distribution and rate. However, this goes 
beyond the scope of the present paper. Here we set con- 
servative detection efficiencies of the data analysis pro- 
cedures. 

The main progresses over the first IGEC analysis 01 
rely on the optimization with respect to amplitude and 
direction of the gw, on the balance between false dis- 
missal probability and false alarm background, on the 
assessment of the statistical coverage of the upper lim- 
its and on the more extended observation time. The 
recent findings from the NAUTILUS and EXPLORER 
2001 data |2(| cannot be compared with the results of 
the analysis presented in this paper, as they relate to 
detectors runs with higher sensitivity. 

We remind once more the reader that we are giving to 
the confidence intervals a frequentist statistical interpre- 
tation: they are determined to the best of our knowledge 
in order to include the actual value of the detected gw 
rate with a relative frequency given on average by the 
chosen coverage. The method and the likclyhood from 
which the confidence intervals are computed is described 
in Section Till El An alternative approach based on the 
Bayesian framework can be followed for the data analy- 
sis. Section IlII El gives all the necessary information, in 
particular the most credible intervals assuming a uniform 
prior. These credible intervals would result similar to the 
ones we presented, but the value for their confidence and 
its interpretation would be different. For instance, the 
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presented confidence intervals with 90% coverage would 
also be the most credible intervals with degree of belief 
of 94%, but different priors can be used as well. 

IGEC amplitude data are in terms of the Fourier com- 
ponent H of the strain amplitude h of the gw burst sig- 
nal. The relation between h and H depend on the spe- 
cific model of signal shape. For instance, for a signal 
consisting of one sinusoidal cycle of 1 ms period h ~ 
(2 • 10 3 Hz) ■ H . For such a source located at the Galac- 
tic Center and emitting isotropically, the estimated mass 
converted in gw would be ~ 0.01M S • (H/1.5 • 10" 21 s) 2 . 

The exchanged IGEC data do not allow the measure- 
ment of the light travel time of gw signals among the 
detector sites. In particular, past IGEC observations 
cannot resolve the gw direction. This capability will 
be achieved when upgraded detectors will demonstrate 
wider sensitivity bandwidths, much greater than 10Hz, 
and will implement signal acquisition and filtering en- 
suring sub-millisecond time resolution|23j. In fact, un- 
der these conditions, we expect that the overall uncer- 
tainties in estimated arrival time of the events will be 
much smaller than the light travel time between detec- 
tor sites. Progress in this respect has been recently 
achieved[2j1 and are expected for the next runs of bar 
detectors |27il28j. 
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APPENDIX A: TIME SERIES AND 
CORRELOGRAMS 



ized delays between points from different series a cross- 
correlogram, and a self-correlogram if t\ = t\. 

This representation has several advantages over a 
straightforward simple histogram of the time delays be- 
tween successive points when it is applied to a Poisson 
point process. A simple time delay histogram describes 
the first order statistics and, for a perfectly homogeneous 
Poisson point process, would be fit by an exponential den- 
sity function. However, this no longer holds true if the 
rate of the point process varies with time. In addition, 
any phase correlation at time lags longer than the average 
time separation is smeared out and not easily identified. 
A self-correlogram, on the other hand, retains much more 
information about the auto-correlation of the time series, 
and is quite insensitive to fluctuations in the rate of the 
point process that occur on timescales longer than the 
range of the correlogram. 

For instance, consider a non-homogeneous Poisson 
point process, with a rate Ai for the first half of the 
time, and a rate A2 for the other half. Its first-order 
delay histogram is the sum of two different exponential 
distributions, while the self correlogram is still flat (be- 
ing the sum of two flat distributions). If on top of the 
random process we add a periodic series with constant 
rate X p <C Ai, A2, then it would be barely discernible in 
the first-order histogram, while in the self-correlogram it 
would appear as a sharp peak at 1/A p . 

The cross-correlation of two time series is related to the 
expected background due to accidental coincidences. If a 
bin of the cross-correlogram has width dt and is centered 
at time lag At, then the counts inside this bin are propor- 
tional to the number of coincidences one would find after 
a time-shift of one series by ±At and with time window 
aperture dt. 

When the cross-correlogram is not flat in a certain 
range of lags, this means that the statistics of the es- 
timated background coincidence counts is not Poisson. 
Deviations from flatness would suggest a correlation of 
the rate of the point processes, either because of a com- 
mon drive acting on more than one detector or because 
of different local signals but with the same periodic char- 
acteristics. A similar feature could also be caused by a 
border effect due to the time fragmentation of the data 
series. 



Given two uniform and uncorrclatcd random point pro- 
cesses -ft, ) . and (to) . „, the set of ordered time differ- 

ences \ A y = t\ — to r between all possible pairs is 
I J t*>t| 

described by a uniform distribution. This can be dis- 
played by projecting all values in {A y } into an his- 
togram. The bin content is the total number of points 
of the original process {t2} igN that fall within a range 
(equal to the bin size) at a specified time lag from each 
of the points in the series {*i} iGN - As the two series 
are independent, this number is a Poisson variable, and 
its average value is bilinear in the number of events in 
{ti} ieN and {t 2 } igN - We call the histogram of general- 



APPENDIX B: TESTING COINCIDENCES AND 
FALSE DISMISSAL 

In order to check if event times and amplitudes in dif- 
ferent detectors are consistent with a common gw signal, 
we use the Bienayme's inequality pflj . 

We recall that for a random variable x with mean 77, 
the absolute value of the residual of x with respect to r\ 
is greater than e with probability P, given by 

P{\x- V \>e}<^lt} (B1) 
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where E { \x — 77 | n } is the n-th absolute central moment 
of x. The Tchebyschcff's inequality is a special case for 
n = 2. This inequality holds true for any statistics of x, 
as long as the moments exist. 

When testing that 77 is the mean value, P is the con- 
servative probability of false dismissal, i.e. of rejecting 
the hypothesis even though it is true. Concerning our 
analysis, we invert this inequality to compute e given P, 
choosing the most convenient order n. e can be conve- 
niently expressed in terms of the standard deviation and 
a non-dimensional multiplier, i.e. e = ka . 

The IGEC exchanged data provide the variance of the 
estimated arrival time of the events as well as the cen- 
tral moments of 2 nd and 4 th order of the amplitude noise 
distribution. To test the consistency of two values mea- 
sured with different detectors we apply the Bienayme's 
inequality to the random variable Xi—Xj , which, in case of 
events generated by the same gw excitation, has a zero 
mean, variance ll$ = cr? = a? + er 2 and 4th moment 
fiyy = /i- 4) + fj,y + 6af(Tj, being of and fxf' the 2 nd and 
4 th order moments of the i-th detector. Then, Eq. IB II 
reduces to P {\xi — Xj\ > kaij} < k~ n fi[^/a"j. and the 
coincidence test \xi — Xj\ < kuij is passed with the re- 
quired maximum false dismissal probability P if we set 

k = \l (fjff I c^) j P , where the order n is used, which 

gives the most stringent check (given that moments of 
order n > 2 are available). 

Regarding the value of the false dismissal P, there is 
an optimal choice to maximize the chances of gw detec- 
tion. The value Pt used for comparison of event times, 
cq. [3] and was chosen in order to balance between the 
conservative probability of detection of a gw coincidence, 
1 — Pt, and the related accidentals, whose number is pro- 
portional to the time window used for the coincidence 

— 1/2 

search and therefore to P T ■ High values for Pt make 
the coincidence search less efficient, since its efficiency 
decreases more rapidly than the related background. On 
the contrary, setting Pt significantly below 5% has the 
drawback of gaining too little in terms of detection prob- 
ability, while increasing significantly the expected false 
alarms and their fluctuations. The choice Pt = 30% 
would maximize precisely the ratio of the (conservative) 
detection efficiency and the corresponding number of ac- 
cidental coincidences. However, what one should really 
care is the uncertainty due to Poisson fluctuations of the 
background rather than the average background itself. 
The contribution of these fluctuations when setting a 
confidence interval depends in a weaker way on Pt (ap- 
proximately it is given by the square root of the average 
background). In the end, this would favor lower values 
for Pt, and the optimal choice is expected to be around 
5 — 10%. These expectations were confirmed by drawing 
the final results for different choices of Pt- 



APPENDIX C: ANALYTIC TABLES OF 
COINCIDENCE COUNTS AND BACKGROUND 
ESTIMATES 

In the following pages the reader can find more de- 
tailed information on the many analyses performed in 
this IGEC search. These trials differ for the directional 
search (none in Tab. [HJ optimization for the Galactic 
Center direction in Tab. IIIip . for the investigated values 
of the gw amplitude threshold H t , and for the configu- 
rations of the network (there where at most 18 disjoint 
choices of detectors, plus the total lines, which synthctize 
the results of all configurations at the same threshold). 
For each trial we report the observation time, the cor- 
responding coincidence counts and the background esti- 
mates. From these data, and assuming that the bursts 
can be modelled as a Poisson point process, we com- 
pute the confidence intervals on the average number of 
gw bursts that possibly occurred within that time span 
(see Section [ill E|) . They are reported in the last two 
columns of the tables for two different values of the min- 
imum coverage, namely 90% and 95%. 

Almost all of the computed confidence intervals cover 
the null result and therefore, at first glance, there is no 
strong evidence for gw detection. In order to be quantita- 
tive in this conclusion, we should undergo the not-so-easy 
task of estimating how many false detections we should 
expect in the tables. To do this, one has to consider the 
exact coverage of the confidence intervals for the spe- 
cific case of no gw present in the data, rather than the 
stated conservative coverage, which is the minimum en- 
sured coverage over any possible number of detected gws. 
Secondly, one must also understand which lines in the ta- 
bles are really independent and which ones are not. 

As already remarked in Section IIV Al the false alarm 
probability is much smaller than the conservative false 
dismissal. Considering the results on the whole, the num- 
ber of detections found are in agreement with the false 
detections, predicted by summing up the expected false 
alarm probability of the hundreds of lines in the tables. 
We must be careful in drawing the conclusions, because 
the trials we are summing on are not completely indepen- 
dent. The following remarks help to get an idea of the 
degree of correlation among different lines in the Tables. 

1. The configurations of detectors at the same search 
threshold H t in each Table arc independent, since 
they refer to mutually disjoint time spans. 

2. Moving toward high H t , as soon as no coincidences 
are found and the observation time saturates to 
100%, the results of each configuration do not de- 
pend anymore on H t (for this reason we simplified 
the Table HJ at high H t ). 

3. The correlation among outcomes found at different 
amplitudes can be anything from zero to one. For 
instance, in Table ITT1 stepping from 3.16 to 3.55 ■ 
10 -20 H z _1 the total coincidence count is the same, 
N c = 8, but only 2 are in common. 
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4. The correlation between the single configuration 
outcomes and the related total line can also vary, 
depending on the relative weight of the configu- 
rations in the sum. In case the total background 
counts are mostly due to a particular configura- 
tion, then the total is quite correlated with it, and 
therefore the total does not to add much informa- 
tion. On the other hand, when a few expected 
total counts are spread in a balanced way among 
all configurations, then each of them may well end 
up showing a coincidence even if its specific back- 
ground is low. These configurations count as many 
independent trials. Therefore, one should not make 
the mistake of selecting just those where a coinci- 
dence was found, because this would lead to a bi- 
ased result. Instead, only the total line for that 
threshold should be considered in this case. 



For all these reasons, the sparse hints of detection that 
can be found here and there in the tables have to be 
criticized from a statistical point of view. In particular, 
we shall discuss two specific cases in more detail. 

The only positive result in a total line of Table ITT1 is at 
H t = 3.98 • 10~ 21 , and seems to confirm the detection 
suggested by the line ex-ni at the same amplitude. How- 
ever, when correctly computing the conditioned probabil- 
ity that the total line is positive when iV c (EX-Nl) = 14 is 
observed, we find it is as high as 54% at 95% coverage. 
Therefore we shall not claim a detection for this total line 
more than we would because of the single ex-ni result, 
and we are already aware of its small significance, con- 
sidering the overall number of expected false detections 
over the single configurations. 

In Tablelinithe only detection are at H t = 3.98 ■ 10~ 21 : 
two detections at different configurations for 0.9 coverage 
and a detection at the total line for both 0.9 and 0.95 
coverage. In order to compute the the probability of get- 
ting a similar result by chance, we devised an empirical 
method [30l | based on the information available from the 
time shifts estimates of the accidental coincidences (see 
Section ITll Df) . We built a thousand independent tables 
of results obtained from the samples of time shifted data 
and we used the resulting statistics of detections as ref- 
erence for the false alarms. The overall number of detec- 
tions confirm the predictions reported in Section llV Af) . 
Moreover, the false alarm probability to show detections 
in at least one total line is 0.55 and 0.33 for 0.9 and 0.95 
coverage respectively. The probability of getting at least 
two detections at the same threshold value is 0.15 at 0.9 
coverage. All probabilities are therefore well in agree- 
ment with no gw detection. 



TABLE II: Results of the time coincidence search analysis for all the 
disjoint configurations of detectors whose observation time has been 
greater than Iday. No directional search has been implemented. The 
coincidence search has been performed with a conservative false dis- 
missal of 5%. The estimated average background Nf, Q$y D is reported 
together with its uncertainty, given either as 5% and 95% percentiles 
or as an 90% upper bound. When followed by the symbol « , the bound 
has to be read as a conservative upper limit (the actual upper bound 
is lower). Where the symbol — appear, there is no available esti- 
mate of the average background, though it docs not affect the related 
confidence intervals because in these cases there were no found coinci- 
dences. For each search threshold, the bottom total line summarizes 
the results, and is obtained by summing over all configurations the ob- 
servation time T ;, s , the number of coincidences 7V C , the background 
Nf, 95%- In the l as t two columns, the confidence intervals on the aver- 
age number of detected gw are reported for 90% and 95% conservative 
coverage. Above the search threshold value 3 ■ 10 — 20 Hz~ 1 the results 
do not change. 
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3 8 





< 1 8 ■ 10 — 4 


— 29 


— 36 




2.3 





< 1.8 • 10 -3 


4 2.9 


4- 3.6 


+ n + nl 

LU LLLL 


558 6 


30 


9Q_l_n o« 


— 13 


— 15 


5.62 


AL-NA 


201.6 





8.4+;^ • 10~ 2 

— 1.5 


4 2.9 


4 3.6 


AL-AU 


85.3 





< 1.8 ■ 10 -3 


4 2.9 


4- 3.6 


AL-EX 


71.5 





~ 1 -1-0 01 Q 
U --0.02 


4 2.9 


4 3.6 


EX-NI 


55.2 


1 


u - ' -0.13 


4 4.7 


4- 5.7 


EX-NA 


36.9 


20 


1 Q 4 + ' 23 
1J -^-0.23 


4 11 


4- 12 


AL-EX-NA 


31.2 





< 1.8 • 10" 4 * 


4 2.9 


4- 3.6 


AL-AU-NA 


28.7 





< 1.8 ■ 10~ 4 * 


4 2.9 


4 3.6 


AL-AU-EX 


26.9 





< 1.8 • 10~ 4 m 


4 2.9 


4 3.6 


AU-EX 


17.0 





1 n 11 'T - ^ 
-'-■"-0.68 lu 


4 2.9 


4 3.6 


AL-AU-NI 


11.2 





< 1.8 • 10~ 4 , 


4 2.9 


4 3.6 


AL-NI 


10.2 





< 1.8 • 10~ 3 » 


4 2.9 


4 3.6 


AL-AU-EX-NA 


9.0 







4 2.9 


4- 3.6 


AL-AU-EX-NI 


8.3 


II 




4 2.9 


4- 3.6 


AU-NI 


6.2 





1 17+O-65 m-2 

1.7I 7 ■ 10 


4 2.9 


4- 3.6 


ATT T7V T\TT 
J\ U HA 11 1 


5 6 





< 1 8 • 10 — 4 


— 29 


— 36 


AT [ ' V T1TT 
AlrUA 111 


5 





< 1 8 • 10 — 4 


— 29 


— 36 




2.4 





< 1.8 • 10 -3 


4 2.9 


4- 3.6 




612 4 


24 


26 + 27 


— 98 


— 12 


6.31 


AL-NA 


198.4 





7.4+J; 4 • 10" 2 


4 2.9 


4 3.6 


AL-EX 


84.5 





r\ 1 +0. 017 
U - J --0.018 


4 2.9 


4 3.6 


AL-AU 


80.4 





< 1.8 ■ 10~ 3 


4 2.9 


4 3.6 


EX-NI 


61.9 


4 


1 + - 13 

°* J -— 0.13 


4 4.9 


4- 5.9 


EX-NA 


45.1 


10 


12.81° 19 


4 6.5 


4- 7.8 


AL-EX-NA 


41.0 





< 1.8 • 10 -4 » 


4 2.9 


4- 3.6 


AL-AU-EX 


29.0 





< 1.8 • 10~ 4 , 


4 2.9 


4 3.6 


AL-AU-NA 


27.2 





< 1.8 • 10~ 4 » 


4 2.9 


4 3.6 


AU-EX 


16.0 





8± 4 ; 4 ■ 10~ 3 


4 2.9 


4 3.6 


AL-NI 


12.8 





< 3.4 ■ 10~ 3 


4 2.9 


4- 3.6 


AL-AU-NI 


11.6 





< 1.8 • 10~ 4 » 


4 2.9 


4- 3.6 


AL-AU-EX-NA 


11.2 







4 2.9 


4- 3.6 


AL-AU-EX-NI 


10.2 







4 2.9 


4 3.6 


AU-EX-NI 


8.6 





< 1.8 • 10 -4 » 


4 2.9 


4 3.6 




6.9 





1 1+ 52 . in- 2 

1 - i -0.56 iU 


— 29 


— 36 


AT [ ' V NT 
rAlj HA .1 1 


6 





< 1 8 • 10 — 4 


— 29 


— 36 




2 3 


q 


< 1 8 • 10 — 3 


— 29 


— 36 


L U L LLL 


653 


l_l 


1 Q + n oil 


— 67 


— 8.1 


7.08 


AL-NA 


193.8 





5.5+\[l ■ 10~ 2 


4 2.9 


4 3.6 


AL-EX 


96.4 





rj 1 +0. 018 
U -- L -0.018 


4 2.9 


4 3.6 


AL-AU 


76.9 





< 1.8 • 10" 3 * 


4 2.9 


4- 3.6 


EX-NI 


65.3 


3 


°- -0.096 


4-5 


4 6 


EX-NA 


49.4 


(i 


n +0 ' 13 

D ' u -0.13 


4 6.5 


4- 7.7 


AL-EX-NA 


48.2 





< 1.8 • 10 -4 » 


4 2.9 


4 3.6 


AL-AU-EX 


32.0 





< 1.8 • 10~ 4 „ 


4 2.9 


4 3.6 


AL-AU-NA 


26.1 





< 1.8 • 10~ 4 , 


4 2.9 


4 3.6 


AU-EX 


16.5 





< 1.8 ■ 10~ 3 


4 2.9 


4- 3.6 


AL-NI 


12.9 





< 3.4 ■ 10 -3 


4 2.9 


4- 3.6 


AL-AU-EX-NA 


12.6 







4 2.9 


4- 3.6 


AL-AU-NI 


11.1 





< 1.8 • 10 4 » 


4 2.9 


4- 3.6 


AL-AU-EX-NI 


10.9 







4 2.9 


4 3.6 


AU-EX-NI 


9.3 





< 1.8 • 10~ 4 „ 


4 2.9 


4 3.6 


AL- EX-NI 


6.8 





< 1.8 • 10 -4 » 


4 2.9 


4- 3.6 


AU-NI 


6.5 





3± 2 - 6 • lO" 3 


4 2.9 


4- 3.6 



continue ^> 
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H-10 21 

(Hz" 1 ) 


config 


T obs 
(days) 


Av 


95% 


confidcnc 

nno/ 
90 /o 


c interval 

90 /o 






2.1 





< 1.8 ■ 10~ 3 


4 2.9 


4 3.6 


LULU,!, 


676 8 


<) 


Q fi + n i a 


— 7.2 


— 85 


7.94 


AL-NA 


189.2 





3 g + 0.96 . 10 -2 


4- 2.9 


4 3.6 


AL-EX 


106.1 





6At\ l ■ 10~ 2 


4- 2.9 


4 3.6 


AL-AU 


74.4 





< 1.8 ■ 10~ 3 » 


4- 2.9 


4 3.6 


BX-NI 


68.2 


(i 


_ ,4_n nsi 

^' -0.081 


0.12 4- 9.3 


4 10 


AL-EX-NA 


54.3 





< 1.8 ■ 10~ 4 * 


4- 2.9 


4 3.6 


EX-NA 


51.6 


2 


n c+0.084 


4 4.4 


4 5.3 


AL-AU-EX 


34.8 





< 1.8 ■ 10~ 4 , 


4 2.9 


4 3.6 


AL-AU-NA 


25.1 





< 1.8 ■ 10~ 4 , 


4 2.9 


4 3.6 


AU-EX 


16.9 





< 1.8 ■ 10~ 3 » 


4 2.9 


4 3.6 


AL-AU-EX-NA 


13.6 







4 2.9 


4 3.6 


AL-NI 


12.7 





< 1.8 ■ 10 -3 » 


4 2.9 


4 3.6 


AL-AU-EX-NI 


11.3 







4 2.9 


4 3.6 


AL-AU-NI 


10.8 





< 1.8 ■ 10~ 4 , 


4 2.9 


4 3.6 


AU-EX-NI 


9.8 





< 1.8 ■ 10~ 4 * 


4 2.9 


4 3.6 




7 1 


o 


< 1 8 ■ 10 — 4 


— 29 


— 36 




6 


q 


< 1.8 ■ 10 — 3 


— 29 


— 36 




2.0 


o 


< 1.8 ■ 10~ 3 * 


4 2.9 


4 3.6 


LtJ LLLL 


694 3 


$ 


5 140 12 


— 92 


— 11 


8.91 


AL-NA 


187.5 





3-llo. 9 4-10 


4 2.9 


4 3.6 


AL-EX 


110.6 





4.9tJ'2 ' 10 ~ 2 


4 2.9 


4 3.6 


AL-AU 


73.8 





< 1.8 ■ 10~ 3 » 


4 2.9 


4 3.6 


EX-NI 


69.5 


2 


_40.067 

' — 0.067 


4 4.8 


4 5.7 


AL-EX-NA 


56.8 





< 1.8 ■ 10~ 4 , 


4 2.9 


4 3.6 


EX-NA 


52.5 


2 


1-6-0.067 


4 4.8 


4 5.8 


AL-AU-EX 


35.8 





< 1.8 ■ 10~ 4 » 


4 2.9 


4 3.6 


AL-AU-NA 


24.8 





< 1.8 ■ 10~ 4 » 


4 2.9 


4 3.6 


AU-EX 


17.1 





4±1. 4 • 10~ 3 


4 2.9 


4 3.6 


AL-AU-EX-NA 


14.1 







4 2.9 


4 3.6 


AL-NI 


12.7 





< 1.8 ■ 10~ 3 » 


4 2.9 


4 3.6 


AL-AU-EX-NI 


11.4 







4 2.9 


4 3.6 


AL-AU-NI 


10.7 





< 1.8 ■ 10~ 4 » 


4 2.9 


4 3.6 


AU-EX-NI 


10.0 





< 1.8 ■ 10~ 4 » 


4 2.9 


4 3.6 




7 2 


o 


< 1 8 ■ 10 — 4 


— 29 


— 36 




5 9 


o 


< 1.8 ■ 10 — 3 


— 29 


— 36 




2 


o 


< 1 8 ■ 10 — 3 


— 29 


— 36 


LtJ LLLL 


702 4 


4 


ij.ixu.UtfO 


— 6 


— 71 


10.00 


AL-NA 


187.2 





z -*-0.S3 lu 


4 2.9 


4 3.6 


AL-EX 


113.0 





4±i j • 10 -2 


4 2.9 


4 3.6 


AL-AU 


73.5 





< 1.8 ■ 10~ 3 , 


4 2.9 


4 3.6 


EX-NI 


70.0 


1 


1 j-i+0.052 
A,u — 0.053 


4 3.9 


4 4.8 


AL-EX-NA 


57.6 





< 1.8 ■ 10~ 4 * 


4 2.9 


4 3.6 


EX-NA 


52.7 


1 


o+O 056 


4 3.9 


4 4.7 


AL-AU-EX 


36.2 





< 1.8 ■ 10~ 4 , 


4 2.9 


4 3.6 


AL-AU-NA 


24.7 





< 1.8 ■ 10~ 4 » 


4 2.9 


4 3.6 


AU-EX 


17.1 





< 1.8 ■ 10 -3 


4 2.9 


4 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4 3.6 


AL-NI 


12.7 





< 1.8 ■ 10~ 3 , 


4 2.9 


4 3.6 


AL-AU-EX-NI 


11.4 







4 2.9 


4 3.6 


AL-AU-NI 


10.7 





< 1.8 ■ 10~ 4 » 


4 2.9 


4 3.6 


AU-EX-NI 


10.0 





< 1.8 ■ 10~ 4 * 


4 2.9 


4 3.6 


AL- EX-NI 


7.2 





< 1.8 ■ 10~ 4 » 


4 2.9 


4 3.6 


AU-NI 


5.9 





< 1.8 ■ 10~ 3 , 


4 2.9 


4 3.6 


AU-NA 


2.0 





< 1.8 ■ 10~ 3 , 


4 2.9 


4 3.6 


total 


706.4 


2 


2.3±0.078 


4 4.5 


4 5.4 


11 ?.? 


AL-NA 


187.1 





1 o + O. 57 , J-.-2 
1 -°-0.61 lv 


4 2.9 


4 3.6 



continue ^> 



H-10 21 
(Hz- 1 ) 


config 


(days) 


AV 


N,, IXv 
u 95% 


confidcnc 
90% 


c interval 
95% 




AL-EX 


113.7 





2.8±°;|^-10 


4 2.9 


4 3.6 


AL-AU 


73.6 





< 1.8 • 10~ 3 , 


4 2.9 


4- 3.6 


EX-NI 


70.1 





„ . -1-0 034 
u ' -0.034 


4 2.9 


4 3.6 


AL-EX-NA 


57.8 





< 1.8 • 10~ 4 * 


4 2.9 


4- 3.6 


EX-NA 


52.9 


1 


-i n - ! - '-' ■ 0^ 1 
A ' u _ 0.051 


4-4 


4 4.8 


AL-AU-EX 


36.3 





< 1.8 • 10" 4 * 


4 2.9 


4- 3.6 


AL-AU-NA 


24.7 





< 1.8 • 10" 4 * 


4 2.9 


4- 3.6 


AU-EX 


17.1 





1.7Z.Q 7 ■ 10 


4 2.9 


4 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4 3.6 


AL-NI 


12.7 





< 1.8 • 10~ 3 , 


4 2.9 


4 3.6 


AL-AU-EX-NI 


11.4 







4 2.9 


4 3.6 


AL-AU-NI 


10.6 





< 1.8 • 10~ 4 » 


4 2.9 


4- 3.6 


AU-EX-NI 


10.1 





< 1.8 • 10~ 4 » 


4 2.9 


4-3.6 


il [ ■ Y MT 


7.2 




< 1 8 • 10 — 4 


— 29 


— 36 




5 8 


o 


< 1 8 • 10 — 3 


— 29 


— 36 




2.0 


o 


< 1.8 • 10~ 3 


4 2.9 


4- 3.6 


LtJ LUjIj 


707 7 


1 


J. . O \JZ u . U D o 


— 38 


— 46 


12.59 


AL-NA 


187.1 





1 1 + - 52 . 1 n -2 
J -- 1 -0.56 1U 


4 2.9 


4- 3.6 


AL-EX 


113.7 





J -- 1 -0.56 1U 


4 2.9 


4 3.6 


AL-AU 


73.6 





< 1.8 • 10" 3 * 


4 2.9 


4 3.6 


EX-NI 


70.2 





0-2-0.026 


4 2.9 


4 3.6 


AL-EX-NA 


57.8 





< 1.8 • 10" 4 * 


4 2.9 


4- 3.6 


EX-NA 


52.9 





„ —-1-0 044 
u - ' -0.045 


4 2.9 


4- 3.6 


AL-AU-EX 


36.4 





< 1.8 • 10~ 4 » 


4 2.9 


4- 3.6 


AL-AU-NA 


24.7 





< 1.8 • 10~ 4 , 


4 2.9 


4 3.6 


AU-EX 


17.2 





-■-■ D -0.6S 1U 


4 2.9 


4 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4 3.6 


AL-NI 


12.7 





< 1.8 • 10" 3 * 


4 2.9 


4- 3.6 


AL-AU-EX-NI 


11.4 





— 


4 2.9 


4- 3.6 


AL-AU-NI 


10.6 





< 1.8 • 10~ 4 * 


4 2.9 


4- 3.6 


AU-EX-NI 


10.1 





< 1.8 • 10~ 4 * 


4 2.9 


4 3.6 


AL-EX.-NI 


7 2 


o 


< 1 8 • 10 — 4 


— 29 


— 36 




5 8 


o 


< 1 8 • 10 — 3 


— 29 


— 36 




2 




< 1 8 - 10 — 3 


— 29 


— 36 




707 8 


q 




— 29 


— 36 


14.12 


AL-NA 


187.1 





3I3 6 • 10~ 3 


4 2.9 


4 3.6 


AL-EX 


113.7 





Sits ■ 10 " 3 


4 2.9 


4 3.6 


AL-AU 


73.6 





< 1.8 • 10" 3 * 


4 2.9 


4 3.6 


EX-NI 


70.2 





„ , _i_n 1 q 

U,J --0.019 


4 2.9 


4- 3.6 


AL-EX-NA 


57.8 





< 1.8 • 10~ 4 * 


4 2.9 


4- 3.6 


EX-NA 


52.9 





u - 4: -0.033 


4 2.9 


4- 3.6 


AL-AU-EX 


36.4 





< 1.8 • 10~ 4 * 


4 2.9 


4 3.6 


AL-AU-NA 


24.7 





< 1.8 • 10" 4 * 


4 2.9 


4 3.6 


AU-EX 


17.2 





< 1.8 • 10" 3 * 


4 2.9 


4 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4- 3.6 


AL-NI 


12.7 





< 1.8 • 10" 3 * 


4 2.9 


4- 3.6 


AL-AU-EX-NI 


11.4 







4 2.9 


4- 3.6 


AL-AU-NI 


10.6 





< 1.8 • 10" 4 * 


4 2.9 


4 3.6 


AU-EX-NI 


10.1 





< 1.8 • 10~ 4 * 


4 2.9 


4 3.6 


AL-EA.-NI 


7 2 


o 


< 1 8 • 10 — 4 


— 29 


— 36 




5 8 




< 1 8 • 10 — 3 


— 29 


— 36 




2 


q 


< 1 8 • 10 — 3 


— 29 


— 36 


— ££t — 


707.9 


() 


0.5 5 ±0.038 


4 2.9 


4-3.6 


15.85 


AL-NA 


187.1 





< 1.8 ■ 10~ 3 


4 2.9 


4 3.6 


AL-EX 


113.7 





6^4.2 ■ 10_3 


4 2.9 


4^ 3.6 


AL-AU 


73.6 





< 1.8 • 10~ 3 » 


4 2.9 


4- 3.6 


EX-NI 


70.2 





6.9±}; 3 • 10 -2 


4 2.9 


4- 3.6 



continue ^> 



H-10 21 

(Hz" 1 ) 


config 


T obs 
(days) 


N c 


" 95% 


confident 
90% 


c interval 
95% 




AL-EX-NA 


57.8 





< 1.8 ■ 10~ 4 » 


4- 2.9 


4 3.6 


EX-NA 


52.9 





r\ Q-f-0.024 
u - -0.024 


4- 2.9 


4 3.6 


AL-AU-EX 


36.4 





< 1.8 ■ 10~ 4 , 


4- 2.9 


4 3.6 


AL-AU-NA 


24.7 





< 1.8 ■ 10~ 4 , 


4- 2.9 


4 3.6 


AU-EX 


17.2 





< 1.8 ■ 10~ 3 » 


4 2.9 


4 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4 3.6 


AL-NI 


12.7 





< 1.8 ■ 10~ 3 , 


4 2.9 


4 3.6 


AL-AU-EX-NI 


11.4 







4 2.9 


4 3.6 


AL-AU-NI 


10.6 





< 1.8 ■ 10~ 4 » 


4 2.9 


4 3.6 


AU-EX-NI 


10.1 





< 1.8 ■ 10~ 4 » 


4 2.9 


4 3.6 


AL-EA.-NI 


7 2 


o 


< 1 8 ■ 10 — 4 


— 29 


— 36 




5 8 


o 


< 1 8 ■ 10 — 3 


— 29 


— 36 




2 


q 


< 1 8 ■ 10 — 3 


— 29 


— 36 


t U t Lib 


707 9 


q 




— 29 


— 36 


17.78 


AL-NA 


187.1 





< 1.8 ■ 10~ 3 , 


4 2.9 


4 3.6 


AL-EX 


113.7 





5± 3 ; 4 • 10- 3 


4 2.9 


4 3.6 


AL-AU 


73.6 





< 1.8 ■ 10~ 3 » 


4 2.9 


4 3.6 


EX-NI 


70.2 





2.lto.7 3 • 10~ 2 


4 2.9 


4 3.6 


AL-EX-NA 


57.8 





< 1.8 ■ 10~ 4 , 


4 2.9 


4 3.6 


EX-NA 


52.9 





r\ i "HO. 017 
U - J --Q.017 


4 2.9 


4 3.6 


AL-AU-EX 


36.4 





< 1.8 ■ 10~ 4 * 


4 2.9 


4 3.6 


AL-AU-NA 


24.7 





< 1.8 ■ 10~ 4 * 


4 2.9 


4 3.6 


AU-EX 


17.2 





< 1.8 ■ 10~ 3 * 


4 2.9 


4 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4 3.6 


AL-NI 


12.7 





< 1.8 ■ 10~ 3 * 


4- 2.9 


4 3.6 


AL-AU-EX-NI 


11.4 





— 


4 2.9 


4 3.6 


AL-AU-NI 


10.6 





< 1.8 ■ 10~ 4 * 


4 2.9 


4 3.6 


AU-EX-NI 


10.1 





< 1.8 ■ 10~ 4 * 


4 2.9 


4 3.6 




7 2 




< 1 8 ■ 10 — 4 


— 29 


— 36 




5 8 


q 


< 1 8 ■ 10 — 3 


— 29 


— 36 




2 


o 


< 1 8 ■ 10 — 3 


— 29 


— 36 


LU LUL 


707 9 


o 




— 29 


— 36 


19.95 


AL-NA 


187.1 





< 1.8 ■ 10~ 3 » 


4 2.9 


4 3.6 


AL-EX 


113.7 





4±l 4 • 10- 3 


4 2.9 


4 3.6 


AL-AU 


73.6 





< 1.8 ■ 10~ 3 , 


4 2.9 


4 3.6 


EX-NI 


70.2 





— -1-4 1 .. ^ 
7_4 b ■ 10 3 


4 2.9 


4 3.6 


AL-EX-NA 


57.8 





< 1.8 ■ 10~ 4 , 


4 2.9 


4 3.6 


EX-NA 


52.9 





4-7±\\\ ■ 10~ 2 


4 2.9 


4 3.6 


AL-AU-EX 


36.4 





< 1.8 ■ 10~ 4 » 


4 2.9 


4 3.6 


AL-AU-NA 


24.7 





< 1.8 ■ 10~ 4 » 


4 2.9 


4 3.6 


AU-EX 


17.2 





< 1.8 ■ 10~ 3 » 


4 2.9 


4 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4 3.6 


AL-NI 


12.7 





< 1.8 ■ 10~ 3 , 


4 2.9 


4 3.6 


AL-AU-EX-NI 


11.4 







4 2.9 


4 3.6 


AL-AU-NI 


10.6 





< 1.8 ■ 10~ 4 » 


4 2.9 


4 3.6 


AU-EX-NI 


10.1 





< 1.8 ■ 10~ 4 » 


4 2.9 


4 3.6 




7 2 


q 


< 1 8 ■ 10 — 4 


— 29 


— 36 




5 8 


o 


< 1 8 ■ 10 — 3 


— 29 


— 36 




2 


o 


< 1 8 ■ 10 — 3 


— 29 


— 36 


LU LUL 


707 9 


o 


7±i 3 ■ 10 — 2 


— 29 


— 36 


22.39 


AL-NA 


187.1 





< 1.8 ■ 10~ 3 » 


4 2.9 


4 3.6 


AL-EX 


113.7 





3+f 6 • 10~ 3 


4 2.9 


4 3.6 


AL-AU 


73.6 





< 1.8 ■ 10~ 3 » 


4 2.9 


4 3.6 


EX-NI 


70.2 





5± 3 ; 4 • 10~ 3 


4 2.9 


4 3.6 


AL-EX-NA 


57.8 





< 1.8 ■ 10~ 4 , 


4 2.9 


4 3.6 


EX-NA 


52.9 





1.8±° ;?£ • 10- 2 


4 2.9 


4 3.6 


AL-AU-EX 


36.4 





< 1.8 ■ 10~ 4 , 


4 2.9 


4 3.6 



continue ^> 
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H10 21 
(Hz- 1 ) 


config 


(days) 


AV 


N,, lf v 

" 95% 


confident 
90% 


c interval 
95% 




AL-AU-NA 


24.7 





< 1.8 ■ 10~ 4 m 


4 2.9 


4 3.6 


AU-EX 


17.2 





< 1.8 • 10~ 3 * 


4 2.9 


4- 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4 3.6 


AL-NI 


12.7 





< 1.8 • 10" 3 * 


4 2.9 


4- 3.6 


AL-AU-EX-NI 


11.4 







4 2.9 


4 3.6 


AL-AU-NI 


10.6 





< 1.8 • 10~ 4 * 


4 2.9 


4- 3.6 


AU-EX-NI 


10.1 





< 1.8 • 10 -4 * 


4 2.9 


4-3.6 




7.2 


o 


< 1.8 • 10~ 4 


4 2.9 


4 3.6 




5 8 


o 


< 1 8 • 10 — 3 


— 29 


— 36 




2 


o 


< 1 8 ■ 10 — 3 * 


— 29 


— 36 


LU LUjIj 


707 9 


o 


1 S-L_n o-} - i n — 2 


— 29 


— 36 


25.12 


AL-NA 


187.1 





< 1.8 ■ 10~ 3 » 


4 2.9 


4 3.6 


AL-EX 


113.7 





3± 2 ' 6 ■ 10~ 3 


4 2.9 


4 3.6 


AL-AU 


73.6 





< 1.8 • 10~ 3 , 


4 2.9 


4 3.6 


EX-NI 


70.2 





4±3.4 ■ 10 ~ 3 


4 2.9 


4 3.6 


AL-EX-NA 


57.8 





< 1.8 • 10~ 4 » 


4 2.9 


4 3.6 


EX-NA 


52.9 





1 1 + u -->^ , i n ^ 
J -- 1 -0.56 1U 


4 2.9 


4- 3.6 


AL-AU-EX 


36.4 





< 1.8 • 10" 4 * 


4 2.9 


4- 3.6 


AL-AU-NA 


24.7 





< 1.8 • 10" 4 * 


4 2.9 


4 3.6 


AU-EX 


17.2 





< 1.8 • 10" 3 * 


4 2.9 


4 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4 3.6 


AL-NI 


12.7 





< 1.8 • 10~ 3 * 


4 2.9 


4- 3.6 


AL-AU-EX-NI 


11.4 





— 


4 2.9 


4- 3.6 


AL-AU-NI 


10.6 





< 1.8 • 10~ 4 * 


4 2.9 


4- 3.6 


AU-EX-NI 


10.1 





< 1.8 • 10~ 4 * 


4 2.9 


4- 3.6 


AT 1? V NI 


7.2 


o 


< 1 8 • 10 — 4 


— 29 


— 36 




5 8 


o 


< 1 8 • 10 — 3 * 


— 29 


— 36 




2.0 


o 


< 1.8 • 10 -3 


4 2.9 


4- 3.6 


LU LUjIj 


707 9 


o 


T_i_ n . i n — 2 

O Zt U . o 1 X U 


— 29 


— 36 


28.18 


AL-NA 


187.1 





< 1.8 • 10~ 3 » 


4 2.9 


4 3.6 


AL-EX 


113.7 





3t' 2 3 6 ■ 10~ 3 


4 2.9 


4 3.6 


AL-AU 


73.6 





< 1.8 • 10~ 3 , 


4 2.9 


4 3.6 


EX-NI 


70.2 





3+ 2 6 • 10~ 3 


4 2.9 


4- 3.6 


AL-EX-NA 


57.8 





< 1.8 • 10~ 4 » 


4 2.9 


4- 3.6 


EX-NA 


52.9 





9± 4 ;! ■ lO" 3 


4 2.9 


4- 3.6 


AL-AU-EX 


36.4 





< 1.8 • 10~ 4 , 


4 2.9 


4 3.6 


AL-AU-NA 


24.7 





< 1.8 • 10~ 4 „ 


4 2.9 


4 3.6 


AU-EX 


17.2 





< 1.8 • 10~ 3 » 


4 2.9 


4 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4- 3.6 


AL-NI 


12.7 





< 1.8 • 10~ 3 » 


4 2.9 


4- 3.6 


AL-AU-EX-NI 


11.4 







4 2.9 


4- 3.6 


AL-AU-NI 


10.6 





< 1.8 • 10~ 4 , 


4 2.9 


4 3.6 


AU-EX-NI 


10.1 





< 1.8 • 10~ 4 * 


4 2.9 


4-3.6 


AT [ ' V MI 
AL LA 111 


7 2 


o 


< 1 8 • 10 — 4 * 


— 29 


— 36 




5 8 


o 


< 1 8 • 10 — 3 


— 29 


— 36 




2 


o 


< 1 8 • 10 — 3 


— 29 


— 36 


L U L U L 


707 9 


q 


2 7 + 75 • 10 — 2 


— 29 


— 36 


31.62 


AL-NA 


187.1 





< 1.8 • 10~ 3 » 


4 2.9 


4 3.6 


AL-EX 


113.7 





< 1.8 ■ 10~ 3 


4 2.9 


4 3.6 


AL-AU 


73.6 





< 1.8 • 10~ 3 » 


4 2.9 


4- 3.6 


EX-NI 


70.2 





< 1.8 ■ 10 -3 


4 2.9 


4- 3.6 


AL-EX-NA 


57.8 





< 1.8 • 10~ 4 , 


4 2.9 


4- 3.6 


EX-NA 


52.9 





r + 3.4 in — 3 
5 -3.8 ■ 10 


4 2.9 


4- 3.6 


AL-AU-EX 


36.4 





< 1.8 • 10~ 4 „ 


4 2.9 


4 3.6 


AL-AU-NA 


24.7 





< 1.8 • 10~ 4 „ 


4 2.9 


4 3.6 


AU-EX 


17.2 





< 1.8 • 10~ 3 » 


4 2.9 


4- 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4- 3.6 



continue ^> 



20 



H-10 21 
(Hz" 1 ) 


config 


T obs 
(days) 


N c 


N,, lf % 


confident 
90% 


c interval 
95% 




AL-NI 


12.7 





< 1.8 ■ 10~ 3 » 


4- 2.9 


4 3.6 


AL-AU-EX-NI 


11.4 







4- 2.9 


4 3.6 


AL-AU-NI 


10.6 





< 1.8 ■ 10~ 4 , 


4- 2.9 


4 3.6 


AU-EX-NI 


10.1 





< 1.8 ■ 10~ 4 » 


4- 2.9 


4 3.6 


AL-EX-NI 


7.2 





< 1.8 ■ 10~ 4 , 


4 2.9 


4 3.6 


AU-NI 


5.8 


II 


< 1.8 ■ 10~ 3 » 


4 2.9 


4 3.6 


AU-NA 


2.0 





< 1.8 ■ 10~ 3 , 


4 2.9 


4 3.6 


total 


707.9 





2.1±0.59 • 10~ 2 


4 2.9 


4 3.6 



TABLE III: Results of the event search analysis as described in the 
previou s tabl e but optimized for the Galactic Center Direction (see 
Section ITTlTl . 



H-10 21 
(Hz- 1 ) 


config 


T obs 

(days) 


N c 


1Nb 95% 


confidence 
90% 


interval 

95% 


2.24 


AU-NA 


1.6 





8 - 1 ii.i ■ 10-2 


4- 2.9 


+ 3.6 


2.51 


AU-NA 


3.1 





q j+0.017 


4 2.9 


4- 3.6 


2.82 


AL-NA 


5.9 


1 


9+° °? 

u - J -0.05 


+ 4 


+ 4.8 


AL-AU 


5.0 


() 


4+0.032 
— 0.032 


4- 2.9 


+ 3.6 


AU-NA 


3.2 


I) 


u - -0.025 


+ 2.9 


+ 3.6 




14.0 


1 


1 .5 + 0.065 


+ 3.8 


4- 4.6 


3.16 


AL-NA 


61.3 


5 


o .- + 0.098 
-0.099 


+ 7 


+ 8.1 


AL-AU 


36.3 


1 


1 ,- + 0.063 
-0.064 


+ 3.8 


4- 4.6 


AL-AU-NA 


4.5 







+ 2.9 


+ 3.6 


AL-EX 


1.7 





n ,-, + 0.024 
U -^-0.024 


4- 2.9 


+ 3.6 


AU-EX 


1.2 





n .-, + 0.024 
Ulii -0.024 


H- 2.9 


+ 3.6 


total 


104.5 


6 


5.5 + 0.12 


+ 6.8 


+ 8 


3.55 


AL-NA 


90.0 





, .+0.062 

^■^-o.oea 


+ 2.9 


+ 3.6 


AL-AU 


54.5 





ri f? + 0.039 
U - D -0.039 


+ 2.9 


+ 3.6 


AL-AU-NA 


9.7 





< 9.1 • 10~ 4 * 


4- 2.9 


+ 3.6 


AL-EX 


5.1 





-i +0 02 

0-llo.02 


H- 2.9 


+ 3.6 


AU-EX 


3.4 





rx 1+0 033 

u * -0.033 


+ 2.9 


4- 3.6 


EX-NI 


2.0 


2 


1 ,- + 0.063 
i,i> — 0.064 


+ 4.9 


4- 5.8 


EX-NA 


1.7 


2 


_i_0 038 
u -°— 0.038 


+ 5.6 


+ 6.5 


AL-AU-EX 


1.6 





< 9.1 ■ 10~ 4 


4- 2.9 


+ 3.6 


total 


167.6 


4 


4.6+o.n 


+ 5.4 


+ 6.5 


3.75 


AL-NA 


97.8 


2 


1 n + 0.053 
' L * u — 0.053 


+ 5.2 


+ 6.1 


AL-AU 


59.1 





n O + 0.028 
U - J -0.029 


+ 2.9 


4- 3.6 


AL-AU-NA 


10.9 





< 9.1 • 10" 4 * 


4- 2.9 


+ 3.6 


AL-EX 


7.3 





b.ltYl -lO" 2 


4 2.9 


+ 3.6 


AU-EX 


4.9 





Ul4: -0.035 


4- 2.9 


+ 3.6 


EX-NI 


3.5 


1 


9 r- + 0.082 
" £, - ,J -0.082 


4 3.6 


+ 4.4 


AL-AU-EX 


2.9 





< 9.1 • 10~ 4 * 


4 2.9 


4- 3.6 


EX-NA 


2.3 





n 7 +0.044 
u - ' -0.045 


4 2.9 


4- 3.6 


AL-AU-EX-NA 


1.8 







4- 2.9 


+ 3.6 


AL-NI 


1.6 





l.7 + _° f ■ lO" 2 


+ 2.9 


+ 3.6 


AU-NA 


1.0 





2.5±°; s 84 • 10- 2 


+ 2.9 


+ 3.6 


total 


192.8 


3 


5.1±0.12 


+ 4.4 


+ 5.4 


3.98 


AL-NA 


105.1 


3 


n „ + 0.046 
u -°-0.046 


0.0061 + 6.8 


4 7.8 


AL-AU 


59.6 





n 9+0.02 

"■^-0.021 


4- 2.9 


+ 3.6 


AL-AU-NA 


12.2 





< 9.1 ■ 10~ 4 


4- 2.9 


+ 3.6 


AL-EX 


10.4 


1 


9.2t\- 6 6 ■ 10~ 2 


4- 4.5 


+ 5.3 


AU-EX 


6.1 


1 


n q+o.03 

u - J -0.03 


+ 4.3 


+ 5.1 


EX-NI 


6.0 





9 Q + 0S9 
z - J -0.089 


+ 8.7 


+ 10 


AL-AU-EX 


4.4 





< 9.1 • 10~ 4 , 


+ 2.9 


+ 3.6 


EX-NA 


3.0 





„ 7 +0.044 
u - ' -0.044 


4- 2.9 


+ 3.6 



continue ^> 



H-10 21 

(Hz^ 1 ) 


config 


(days) 




u 95% 


confidence 
90% 


interval 

95% 




AL-AU-EX-NA 


2.8 







+ 2.9 


+ 3.6 


AL-NI 


2.3 


1 


o p + 0.95 1 n — 2 
— 0.99 


0.014 4 4.6 


+ 5.4 


4 T AIT MT 
AL-rtU 1^ 1 


1 9 


Q 




— 29 


— 36 


AL-AU-EX-NI 


1.4 


o 




+ 2.9 


4- 3.6 




1 i 


o 


9 ,+0.73 ,„-2 
z,i -0.77 lu 


— 29 


— 36 




215 8 


12 




1 5 — 15 


88 — 16 


4.47 


AL-NA 


116.9 





n e + o ■ 036 
u - lJ -0.036 


4- 2.9 


+ 3.6 


AL-AU 


60.3 





S.6t\'l ■ 10~ 2 


4 2.9 


+ 3.6 


AL-EX 


18.0 





8.5+ 1 ! ■ 10~ 2 


+ 2.9 


4- 3.6 


AL-AU-NA 


14.4 





< 9.1 • 10~ 4 » 


4 2.9 


+ 3.6 


EX-NI 


13.1 


5 


^■ 1 -0.11 


+ 6.6 


4- 7.8 


AL-AU-EX 


7.5 





< 9.1 • 10 -4 * 


+ 2.9 


+ 3.6 


AU-EX 


7.5 





ri i-j + . 02 
U - Z -0.02 


4 2.9 


+ 3.6 


EX-NA 


5.6 


4 


°- U -0.09 


4^ 6.3 


+ 7.4 


AL-AU-EX-NA 


5.1 







+ 2.9 


+ 3.6 


AL-AU-EX-NI 


4.4 







4 2.9 


+ 3.6 


AL-AU-NI 


3.5 







+ 2.9 


+ 3.6 


AL-NI 


3.2 





n Q+0.57 i n — 2 

1.3T 61 ■ 10 


+ 2.9 


+ 3.6 


AL-EX-N A 


2.3 


o 




+ 2.9 


+ 3.6 




1.7 


Q 




— 29 


— 36 




1 2 


o 


3 7+0.98 _ lQ-2 


— 29 


— 36 


total 


264.7 


9 


8+0 15 


+ 8.1 


+ 9.4 


5.01 


AL-NA 


125.3 





fl 4+0.034 
u - -0.035 


+ 2.9 


+ 3.6 


AL-AU 


59.1 


1 


„ .4-1 S ^ „ o 

6.41^3 ■ 10 2 


4- 4.5 


+ 5.4 


AL-EX 


27.5 





7.9tJ- 4 ■ 10" 2 


4 2.9 


+ 3.6 


EX-NI 


21.2 


2 


^' — 0.11 


044 


+ 4.9 


AL-AU-NA 


13.2 







4- 2.9 


+ 3.6 


EX-NA 


9.8 


10 


— -4-fl 15 

7-8lo. 15 


4- 9.2 


4 11 


AL-AU-EX 


9.7 







4- 2.9 


+ 3.6 


AU-EX 


8.6 





^■2-0.023 


4 2.9 


+ 3.6 


AL-AU-EX-NA 


7.3 







4 2.9 


+ 3.6 


AL-AU-EX-NI 


6.3 







4 2.9 


+ 3.6 


AL-EX-NA 


4.7 





— 


4- 2.9 


4- 3.6 


AL-AU-NI 


3.9 







4- 2.9 


+ 3.6 


AL-NI 


3.0 





6 -4.2 ■ 10 ~ 


+ 2.9 


+ 3.6 




2 5 


I 


7 Q+ 1 - 4 . in- 2 


— 45 


— 5.3 


a T T [ ' V MT 


1 6 


o 




— 2.9 


— 36 




1 3 


o 


i 8 +0 - 67 ■ in" 2 


— 2.9 


— 36 


LO Lai 


306 8 


l_l 


13 + 19 


— 96 


— 11 


5.62 


AL-NA 


130.8 





0-2^0.026 


4- 2.9 


4- 3.6 


AL-AU 


57.3 





3.9+ 1 • 10~ 2 


4 2.9 


+ 3.6 


AL-EX 


37.5 





ri 1+0.017 
U - A — 0.018 


4 2.9 


+ 3.6 


EX-NI 


29.4 


(i 




4- 7.5 


+ 8.7 


EX-NA 


15.4 


12 


10.5±° -\l 


4- 9.3 


+ 11 


AL-AU-NA 


14.3 







4- 2.9 


4- 3.6 


AL-AU-EX 


12.7 







4- 2.9 


+ 3.6 


AU-EX 


9.6 





^■■^-0.022 


4 2.9 


+ 3.6 


AL-AU-EX-NA 


9.0 







+ 2.9 


+ 3.6 


AL-AU-EX-NI 


8.3 







4 2.9 


+ 3.6 


AL-EX-NA 


8.3 







4- 2.9 


4- 3.6 


AL-NI 


5.1 





4±| 4 ■ 10 -3 


4 2.9 


+ 3.6 


AL-AU-NI 


4.5 







4- 2.9 


+ 3.6 


AU-NI 


3.1 





2.6I ; 86 ■ 10 2 


4- 2.9 


+ 3.6 


AU-EX-NI 


2.4 







4 2.9 


+ 3.6 


AL-EX-NI 


1.5 







4 2.9 


+ 3.6 


AU-NA 


1.4 





8±tt ■ lO" 3 


4- 2.9 


+ 3.6 


total 


352.9 


IS 


16 + 0.2 


4 11 


+ 13 



continue ^> 



21 



H-10 21 
(Hz- 1 ) 


con fig 


T bs 

(days) 




V, 5% 
iNb 95% 


confidence 
9U /o 


interval 

90/0 


6.31 


AL-NA 


134.7 





_I_q 026 
-0.026 


4 2.9 


4- 3.6 


AL-AU 


55.6 





Lfl Qfi 9 

3.6l"- J ' 10 


4 2.9 


4- 3.6 


AL-EX 


48.0 





ri -| + . 01 8 

U - A — 0.018 


4 2.9 


4- 3.6 


EX-NI 


34.7 


1 


2-9-0.088 


4- 3.5 


4- 4.3 


EX-NA 


20.9 


8 


I0.5t2;£ 


4- 6.1 


4- 7.3 


AL-AU-NA 


15.2 





~~ 


4- 2.9 


4- 3.6 


AL-AU-EX 


15.1 







4- 2.9 


4- 3.6 


AL-EX-NA 


12.4 





~~ 


4 2.9 


4- 3.6 


AL-AU-EX-NA 


11.2 







4- 2.9 


4- 3.6 


AU-EX 


10.2 





ri -| + . 01 8 

U - A — 0.018 


4 2.9 


4- 3.6 


AL-AU-EX-NI 


10.2 





— 


4- 2.9 


4- 3.6 


AL-NI 


6.7 





< 1.8 ■ 10~ 3 


4- 2.9 


4- 3.6 


AL-AU-NI 


5.3 







4 2.9 


4- 3.6 


AU-NI 


3.6 







4- 2.9 


4- 3.6 


ATT T7V MT 
1\ U HA 11 1 


3 3 


o 




— 29 


— 36 


AL-EX-NI 


2 3 


o 




— 29 


— 36 




1 5 




< 1 8 ■ 10 — 3 


— 29 


— 36 


hU lu,l 


392 8 


() 




— 56 


— 68 


7.08 


AL-NA 


137.5 





_[_q 026 
— 0.027 


4 2.9 


4- 3.6 


AL-EX 


57.8 





ri i ~l"0. 019 
U '- L — 0.02 


4 2.9 


4- 3.6 


AL-AU 


54.8 





g ■ 10 -3 


4- 2.9 


4- 3.6 


EX-NI 


39.5 


:s 


_ _ _i_n dr i 
z -°-0.085 


4- 5.4 


4- 6.4 


EX-NA 


26.5 


7 


io.3±S:l7 


4- 5.5 


4- 6.7 


AL-AU-EX 


17.7 







4- 2.9 


4- 3.6 


AL-EX-NA 


16.6 





— 


4 2.9 


4- 3.6 


AL-AU-NA 


15.7 







4- 2.9 


4- 3.6 


AL-AU-EX-NA 


12.6 







4- 2.9 


4- 3.6 


AL-AU-EX-NI 


10.9 





~~ 


4- 2.9 


4- 3.6 


AU-EX 


10.9 





6.5t\;l ■ 10 ~ 2 


4- 2.9 


4- 3.6 


AL-NI 


7.5 





4+| 4 ■ 10~ 3 


4- 2.9 


4- 3.6 


AL-AU-NI 


5.8 







4 2.9 


4- 3.6 


AU-EX-NI 


4.4 







4- 2.9 


4-3.6 




3 5 


o 


9+ 4 7 ■ 10~ 3 


— 29 


— 36 


AT [ ' Y MT 


3 2 


o 




— 29 


— 36 




1 4 


() 


< 1.8 ■ 10™ 3 


— 29 


— 36 


LU L It t 


428 1 


10 


43 1 q 


— 63 


— 7.6 


7.94 


AL-NA 


139.1 





0-2-0.023 


4 2.9 


4- 3.6 


AL-EX 


67.2 





9.6^11 ■ 10~ 2 


4- 2.9 


4- 3.6 


AL-AU 


54.3 





z -0.75 lu 


4- 2.9 


4- 3.6 


EX-NI 


43.4 


4 


2 '3-0.08 


4- 6.7 


4- 7.8 


EX-NA 


30.8 


9 


+ 0. 13 
— 0.13 


4- 9.5 


4-11 


AL-EX-NA 


20.9 







4 2.9 


4- 3.6 


AL-AU-EX 


20.0 







4 2.9 


4- 3.6 


AL-AU-NA 


16.1 







4- 2.9 


4- 3.6 


AL-AU-EX-NA 


13.6 







4- 2.9 


4- 3.6 


AU-EX 


11.5 





6.3±i 3 ■ 10~ 2 


4- 2.9 


4- 3.6 


AL-AU-EX-NI 


11.3 







4 2.9 


4- 3.6 


AL-NI 


8.0 





< 1.8 ■ 10~ 3 


4 2.9 


4- 3.6 


AL-AU-NI 


6.1 





~~ 


4- 2.9 


4- 3.6 


AU-EX-NI 


5.1 







4- 2.9 


4- 3.6 




3 9 


o 




— 29 


— 36 




3 5 


o 


< 1.8 ■ 10 3 


— 29 


— 36 


AU-NA 


1.4 





5+l i ■ 10" 3 

— 3.8 


4- 2.9 


4- 3.6 


total 


458.2 


13 


8.9±o.i6 


4- 12 


4- 13 


8.91 


AL-NA 


141.2 





n .-, + 0.025 
"■^-0.025 


4- 2.9 


4- 3.6 


AL-EX 


75.2 





n q +0.017 
u- — 0.018 


4- 2.9 


4- 3.6 


AL-AU 


55.1 





gj: 4 -.; ■ lO" 3 


4 2.9 


4- 3.6 



continue ^> 



H-10 21 

(Hz^ 1 ) 


config 


T„bs 

(days) 


N 


u 95% 


confidence 
90% 


interval 

95% 




EX-NI 


46.7 


2 


■7+0 068 
I- ' -0.068 


4- 4.8 


4- 5.7 


EX-NA 


33.5 


5 




4- 6.8 


4-8 


AL-EX-NA 


24.9 





~~ 


4 2.9 


4- 3.6 


AL-AU-EX 


21.9 







4- 2.9 


4- 3.6 


AL-AU-NA 


16.5 







4 2.9 


4- 3.6 


AL-AU-EX-NA 


14.1 





_ 


4- 2.9 


4- 3.6 


AU-EX 


12.0 





4.5± i ■ 10~ 2 


4- 2.9 


4- 3.6 


AL-AU-EX-NI 


11.4 







4 2.9 


4- 3.6 


AL-NI 


8.4 







4 2.9 


4- 3.6 


AL-AU-NI 


6.4 





— 


4 2.9 


4- 3.6 


AU-EX-NI 


5.7 







4- 2.9 


4-3.6 


AT [ ' Y MT 

rtij HA 1\ I 


4.4 


o 




— 29 


— 36 




3 7 


o 


+4.4 in -3 
°-4.8 1U 


— 29 


— 36 




1 4 


o 


< 3 4 • 10 — 3 


— 29 


— 36 


tntn 1 

(. U L Ui I 


485 


7 




— 75 


— 88 


10.00 


AL-NA 


144.0 





r\ -I -)-0 . 1 S 
U - 1 -0.018 


4 2.9 


4- 3.6 


AL-EX 


81.2 





ri -I T 1 - 1 . U -LU 

U - A _ 0.017 


4- 2.9 


4- 3.6 


AL-AU 


56.0 





1 4+0-59 , |-,-2 


4- 2.9 


4- 3.6 


EX-NI 


49.3 


1 


, __i_o nfi i 

LO - 0.064 


4- 3.8 


4- 4.6 


EX-NA 


35.6 


2 


r> n + " 089 
^•^-0.09 


4- 4.3 


4- 5.2 


AL-EX-NA 


28.0 







4- 2.9 


4- 3.6 


AL-AU-EX 


23.2 





~~ 


4- 2.9 


4- 3.6 


AL-AU-NA 


16.9 





~~ 


4- 2.9 


4- 3.6 


AL-AU-EX-NA 


14.2 







4- 2.9 


4- 3.6 


AU-EX 


12.6 





5.5+ 1 ! ■ 10 -2 


4- 2.9 


4- 3.6 


AL-AU-EX-NI 


11.4 







4- 2.9 


4- 3.6 


AL-NI 


8.8 





4^3.4 ■ lO"' 


4 2.9 


4- 3.6 


AL-AU-NI 


6.6 







4- 2.9 


4- 3.6 


ALI-EX-NI 


6.1 







4- 2.9 


4- 3.6 


AT [ ' Y MT 

rtij £..Y 1\ I 


4.9 


o 




— 29 


— 36 




3 8 


o 


. + 3.4 1n -3 
°-3.8 1U 


— 29 


— 3.6 




1 5 


o 


< 1 8 ■ 10 — 3 


— 29 


— 36 


6 U L It I 


506 7 


3 




— 45 


— 5.5 


11.22 


AL-NA 


146.8 





i 1 0-019 
u '-0.02 


4- 2.9 


4- 3.6 


AL-EX 


85.0 





5.2tJ; 2 ■ 10~ 2 


4- 2.9 


4- 3.6 


AL-AU 


57.2 





_ ^-l-0fi7 9 

l.sto.n ■ 10 


4 2.9 


4- 3.6 


EX-NI 


51.4 





1 q ~\~ . 06 
1 - J -0.06 


4 2.9 


4- 3.6 


EX-NA 


37.1 


3 


cy /?-|-0.084 

z - u -0.084 


4 5.4 


4- 6.4 


AL-EX-NA 


30.6 





_ 


4- 2.9 


4- 3.6 


AL-AU-EX 


24.0 







4- 2.9 


4- 3.6 


AL-AU-NA 


17.3 







4- 2.9 


4- 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4- 3.6 


AU-EX 


12.9 





xn QQ o 
3 g +0.9J _ 1Q 2 


4 2.9 


4- 3.6 


AL-AU-EX-NI 


11.4 







4 2.9 


4- 3.6 


AL-NI 


9.2 





3t 2 - 6 ■ 10 -3 


4- 2.9 


4- 3.6 


AL-AU-NI 


6.8 





_ 


4- 2.9 


4- 3.6 


ALI-EX-NI 


6.5 







4- 2.9 


4-3.6 


AT [ ' Y MT 
nL 111 


5 2 


o 




— 2.9 


— 36 




3 9 


Q 


a + 3 .in -3 

^-3.4 1U 


— 2.9 


— 36 




1 5 


Q 


< 1 8 ■ 10 — 3 


— 29 


— 36 


LU LU.L 


523 3 


3 


4 2±o n 


— 47 


— 57 


12.59 


AL-NA 


149.3 





rt rj + . 02 1 
u -' 1 -0.022 


4 2.9 


4- 3.6 


AL-EX 


87.7 





n hr+ 1-6 i n — 2 

9-7 i e ■ 10 


4 2.9 


4- 3.6 


AL-AU 


58.3 





■ lO- 2 


4 2.9 


4- 3.6 


EX-NI 


53.0 


3 


,+0-055 
1 - i -0.055 


4- 6.5 


4- 7.5 


EX-NA 


38.3 


1 


, jj+0.07 
^•"-O.OTl 


4- 3.7 


4-4.5 


AL-EX-NA 


32.7 







4- 2.9 


4- 3.6 



continue ^> 



22 



H-10 21 

(Hz" 1 ) 


con fig 


T bs 

(days) 




N 6 %X% 

u 95% 


confidence 
90% 


interval 

95% 




AL-AU-EX 


24.6 







4- 2.9 


4- 3.6 


AL-AU-NA 


17.5 







4- 2.9 


4- 3.6 


AL-AU-EX-NA 


14.2 







4- 2.9 


4- 3.6 


AU-EX 


13.2 





o o-i-u.oo i n 

z -°-0.89 1U 


4- 2.9 


4- 3.6 


AL-AU-EX-NI 


11.4 







4 2.9 


4- 3.6 


AL-NI 


9.5 





< 1.8 ■ 10~ 3 


4- 2.9 


4- 3.6 


AL-AU-NI 


7.0 





— 


4- 2.9 


4- 3.6 


AU-EX-NI 


6.7 







4- 2.9 


4- 3.6 


AL-EA.-NI 


5 3 


o 




— 29 


— 36 




4.1 


o 


< 1.8 ■ 10 — 3 


— 29 


— 36 




1 6 


o 


< 1.8 ■ 10 — 3 


— 2.9 


— 36 




536 4 


4 




— 61 


— 7.2 


14.12 


AL-NA 


151.6 





rj i -J-Q .017 
u *— 0.017 


4- 2.9 


4- 3.6 


AL-EX 


89.7 





7 - 5 -i 1 ■ 10 ~ 2 


4- 2.9 


4- 3.6 


AL-AU 


59.3 





■ 10 ~ 3 


4 2.9 


4- 3.6 


EX-NI 


54.4 


2 


n q+0.048 
u - y -0.048 


4- 5.3 


4- 6.3 


EX-NA 


39.4 


1 


_[_Q 007 

1- ' -0.068 


4- 3.7 


4- 4.5 


AL-EX-NA 


34.4 







4- 2.9 


4- 3.6 


AL-AU-EX 


25.2 





~~ 


4- 2.9 


4- 3.6 


AL-AU-NA 


17.8 





~~ 


4- 2.9 


4- 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4- 3.6 


AU-EX 


13.5 





4.3+JJ ■ 10~ 2 


4- 2.9 


4- 3.6 


AL-AU-EX-NI 


11.4 







4- 2.9 


4- 3.6 


AL-NI 


9.8 





< 1.8 ■ 10~ 3 


4- 2.9 


4- 3.6 


AL-AU-NI 


7.2 







4 2.9 


4- 3.6 


AU-EX-NI 


6.9 







4 2.9 


4- 3.6 


AL-EA.-NI 


5 5 


o 




— 29 


— 36 




4.2 


o 


3 + 2.6 _ 1Q -3 


4- 2.9 


4- 3.6 




1 6 


() 


7+*-l . 1(1-3 
'-4.5 iU 


— 29 


— 36 




547 5 


3 




— 53 


— 63 


15.85 


AL-NA 


153.7 


1 


i-i -I -[-0 . 01 8 
— 0.01 9 


4- 4.5 


4- 5.3 


AL-EX 


91.4 





4.8t}; 2 ■ 10 ~ 2 


4 2.9 


4- 3.6 


AL-AU 


60.2 





6^42 ■ 10 ~ 3 


4- 2.9 


4- 3.6 


EX-NI 


55.9 





n q + 0.049 
u ' 3 — 0.049 


4- 2.9 


4- 3.6 


EX-NA 


40.3 


1 


1-3 — 0.06 


4- 3.8 


4- 4.7 


AL-EX-NA 


35.7 





_ 


4- 2.9 


4- 3.6 


AL-AU-EX 


25.7 







4 2.9 


4- 3.6 


AL-AU-NA 


18.0 







4- 2.9 


4- 3.6 


AL-AU-EX-NA 


14.2 







4- 2.9 


4- 3.6 


AU-EX 


13.7 





3-2io.95 ' 10_2 


4- 2.9 


4- 3.6 


AL-AU-EX-NI 


11.4 







4- 2.9 


4- 3.6 


AL-NI 


10.1 





< 1.8 ■ 10" 3 


4 2.9 


4- 3.6 


AL-AU-NI 


7.4 





— 


4- 2.9 


4- 3.6 


AU-EX-NI 


7.1 







4- 2.9 


4-3.6 




5 6 


() 




— 29 


— 36 




4.3 




3 + 2.6 _ 1Q -3 


— 29 


— 36 




1 6 


o 


< 1 8 • 10 — 3 * 


— 2.9 


— 36 


LU L It t 


557 5 


9 




— 45 


— 54 


17.78 


AL-NA 


155.7 





i-i | ~r0 .018 
u ■ 1 - . 1 9 


4 2.9 


4- 3.6 


AL-EX 


93.0 





4 " 3 — l" 1 ' 10 ~ 2 


4- 2.9 


4- 3.6 


AL-AU 


61.0 





< 1.8 ■ 10" 3 


4- 2.9 


4- 3.6 


EX-NI 


57.0 





fl Q + 0-047 
u -°-0.047 


4- 2.9 


4- 3.6 


EX-NA 


41.2 


1 


-, +0.053 
i --0.054 


4- 3.9 


4- 4.8 


AL-EX-NA 


36.7 







4 2.9 


4- 3.6 


AL-AU-EX 


26.1 







4- 2.9 


4- 3.6 


AL-AU-NA 


18.3 







4- 2.9 


4- 3.6 


AL-AU-EX-NA 


14.2 







4- 2.9 


4- 3.6 



continue ^> 



H-10 21 

(Hz^ 1 ) 


config 


Tab, 

(days) 




u 95% 


confidence 
90% 


interval 

95% 




AU-EX 


13.8 





9 R+ U tt - .in z 


4- 2.9 


4- 3.6 


AL-AU-EX-NI 


11.4 







4- 2.9 


4- 3.6 


AL-NI 


10.4 





< 1.8 • 10" 3 


4 2.9 


4- 3.6 


AL-AU-NI 


7.6 







4- 2.9 


4- 3.6 


AU-EX-NI 


7.2 







4 2.9 


4-3.6 


AT [ ' V MT 


5.7 


o 




— 29 


— 36 




4.5 


o 


< 3 4 • 10 — 3 


— 29 


— 36 




1.7 


o 


< 3.4 • 10~ 3 


4- 2.9 


4- 3.6 


(. U L It I 


566 6 


I 


9 1 _i_n n-7f. 


— 36 


— 4.5 


19.95 


AL-NA 


157.6 


1 


n i + 017 
u ' x — 0.017 


4-4.5 


4- 5.3 


AL-EX 


94.4 





4.2^J 1 ■ 10~ 2 


4- 2.9 


4- 3.6 


AL-AU 


61.8 





5iii ■ lO" 3 


4- 2.9 


4- 3.6 


EX-NI 


58.1 


1 


fl Q~l - Q.045 
u -°-0.046 


4- 4.1 


4- 4.9 


EX-NA 


41.9 





1 n ~t~ . 5 2 
-0.052 


4- 2.9 


4- 3.6 


AL-EX-NA 


37.5 







4- 2.9 


4- 3.6 


AL-AU-EX 


26.5 







4 2.9 


4- 3.6 


AL-AU-NA 


18.5 







4- 2.9 


4- 3.6 


AL-AU-EX-NA 


14.2 







4- 2.9 


4- 3.6 


AU-EX 


14.0 





1 - ± -0.56 ±U 


4- 2.9 


4- 3.6 


AL-AU-EX-NI 


11.4 





_ 


4 2.9 


4- 3.6 


AL-NI 


10.7 







4 2.9 


4 3.6 


AL-AU-NI 


7.7 







4 2.9 


4- 3.6 


ALI-EX-NI 


7.3 







4- 2.9 


4- 3.6 


AT [ ' V MT 


5 8 


o 




— 29 


— 36 




4.6 


o 


^-3.4 ±U 


4- 2.9 


4- 3.6 




1.7 


o 


< 1 8 • 10 — 3 


— 29 


— 36 


6 U L It t 


574.8 


2 




— 4.7 


— 56 


22.39 


AL-NA 


159.3 





4.8+J 1 ■ 10 2 


4- 2.9 


4- 3.6 


AL-EX 


95.7 





7.111;* ■ 


4- 2.9 


4- 3.6 


AL-AU 


62.5 





< 3.4 • 10 -3 


4 2.9 


4- 3.6 


EX-NI 


59.0 


2 


fi+ 04 
U - D -0.041 


4- 5.5 


4- 6.5 


EX-NA 


42.7 


1 


1 j-j-|-0 . 052 
-L-U— 0.053 


4 3.9 


4-4.8 


AL-EX-NA 


38.3 





_ 


4- 2.9 


4- 3.6 


AL-AU-EX 


26.9 







4 2.9 


4- 3.6 


AL-AU-NA 


18.7 







4- 2.9 


4- 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4- 3.6 


AU-EX 


14.2 







4 2.9 


4- 3.6 


AL-AU-EX-NI 


11.4 







4 2.9 


4 3.6 


AL-NI 


11.1 





< 1.8 • 10~ 3 


4- 2.9 


4- 3.6 


AL-AU-NI 


7.8 







4- 2.9 


4- 3.6 


ALI-EX-NI 


7.4 







4- 2.9 


4-3.6 


AT [ ' V MT 

ftL INI 


5 9 


Q 




— 2.9 


— 36 




4.7 


o 


q+4.7 1() -3 
— 5.1 1U 


— 2.9 


— 36 




1.7 


Q 


< 1 8 • 10 — 3 * 


— 2.9 


— 36 


LU LU.L 


582 6 




1 8 + 07 


— 59 


— 69 


25.12 


AL-NA 


161.0 





7.8±\ l ■ 10~ 2 


4- 2.9 


4- 3.6 


AL-EX 


96.9 





7. 1 + 1* ■ 10" 2 


4 2.9 


4- 3.6 


AL-AU 


63.2 





6 -4.2 ■ 10_3 


4 2.9 


4 3.6 


EX-NI 


59.7 


1 


r\ [- + 0.035 
UltJ - 0.035 


4 4.2 


4- 5.1 


EX-NA 


43.3 





1 j-,-j-0 .052 
-"-■ u -0.053 


4- 2.9 


4- 3.6 


AL-EX-NA 


39.0 







4- 2.9 


4- 3.6 


AL-AU-EX 


27.2 







4- 2.9 


4- 3.6 


AL-AU-NA 


18.9 







4- 2.9 


4- 3.6 


AU-EX 


14.4 





n "7+0.83 — 2 
z - ' -0.87 ' lu 


4 2.9 


4- 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4- 3.6 


AL-AU-EX-NI 


11.4 







4- 2.9 


4- 3.6 


AL-NI 


11.4 





< 1.8 • 10 -3 , 


4- 2.9 


4- 3.6 



continue ^> 



H-10 21 

(Hz" 1 ) 


con fig 


T bs 

(days) 




N 6 %X% 


confidence 
90% 


interval 

95% 




AL-AU-NI 


7.9 







4 2.9 


4 3.6 


AU-EX-NI 


7.5 







4 2.9 


4 3.6 


AL-EX.-NI 


6 


o 




— 29 


— 36 




4.8 


o 


^-3.4 1U 


-f- 2.9 


4 3.6 




1 7 


o 


< 1 8 • 10 — 3 


— 29 


— 36 




589 7 


1 


1 7-l-n nR7 

J- . 1 IjZ U.UD i 


— 37 


— 46 


28.18 


AL-NA 


162.6 





r\ i ~HQ -017 
u - ± -0.018 


4 2.9 


4 3.6 


AL-EX 


98.0 





6.6+1 1 ' 10 


4 2.9 


4 3.6 


AL-AU 


63.8 





3+3' 6 • 10~ 3 


4 2.9 


4 3.6 


EX-NI 


60.3 


1 


r\ q -|-0 . 027 
u -°— 0.028 


4 4.3 


4 5.2 


EX-NA 


43.9 





„ „_i_n nfi 


4 2.9 


4 3.6 


AL-EX-NA 


39.6 





~~ 


4 2.9 


4 3.6 


AL-AU-EX 


27.5 





~~ 


4 2.9 


4 3.6 


AL-AU-NA 


19.1 





~~ 


4 2.9 


4 3.6 


AU-EX 


14.6 







4 2.9 


4 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4 3.6 


AL-NI 


11.6 





< 3.4 ■ 10~ 3 


4 2.9 


4 3.6 


AL-AU-EX-NI 


11.4 







4 2.9 


4 3.6 


AL-AU-NI 


7.9 







4 2.9 


4 3.6 


AU-EX-NI 


7.6 







4- 2.9 


4 3.6 


AL-E.X.-NI 


6 


o 




— 29 


— 36 




4.9 


o 


< 3.4 ■ 10~ 3 


4 2.9 


4 3.6 




1.7 


() 


< 1 8 • 10 — 3 


— 29 


— 36 




596 1 


1 


1 A. -J- n n^i 


— 38 


— 46 


31.62 


AL-NA 


164.1 





5.6±}- 2 ■ 10~ 2 


4 2.9 


4 3.6 


AL-EX 


99.0 





7.4tJ * ■ 10~ 2 


4 2.9 


4 3.6 


AL-AU 


64.5 





3+g- 6 • 10~ 3 


4 2.9 


4 3.6 


EX-NI 


60.9 


1 


„ ^ + 024 
"■^-0.024 


4 4.4 


4 5.2 


EX-NA 


44.4 


1 


n o~l"^-0-i7 

u -° -0.047 


4 4 


4 4.9 


AL-EX-NA 


40.2 





_ 


4 2.9 


4 3.6 


AL-AU-EX 


27.8 







4 2.9 


4 3.6 


AL-AU-NA 


19.2 







4 2.9 


4 3.6 


AU-EX 


14.7 





2.11° „ • 10 2 


4 2.9 


4 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4 3.6 


AL-NI 


11.7 





< 1.8 ■ io -3 


4 2.9 


4 3.6 


AL-AU-EX-NI 


11.4 







4- 2.9 


4 3.6 


AL-AU-NI 


8.0 





~~ 


4 2.9 


4 3.6 


AU-EX-NI 


7.7 







4 2.9 


4 3.6 




6 1 






— 29 


— 36 




5 


o 


a+ 3 .in -3 

4 -3.4 iU 


— 29 


— 36 




1 7 


o 


4 +3 -in -3 


— 29 


— 36 


tntn 1 

LU L It t 


602 1 


9 


± . ^ It U . U <J 1 


— 51 


— 6 


35.48 


AL-NA 


165.6 





(-1 1 ~l~0 . 01 8 
u ■ 1 — .0 19 


4 2.9 


4 3.6 


AL-EX 


99.9 





4.6tJ J ■ 10" 2 


4 2.9 


4 3.6 


AL-AU 


65.1 





9±* : ; ■ lO" 3 


4 2.9 


4 3.6 


EX-NI 


61.5 





*-i .-,+0.022 

U - Z -0.022 


4 2.9 


4 3.6 


EX-NA 


45.0 





_ __i_n n 1 i"i 

U- ° -0.046 


4 2.9 


4 3.6 


AL-EX-NA 


40.8 







4 2.9 


4 3.6 


AL-AU-EX 


28.1 







4 2.9 


4 3.6 


AL-AU-NA 


19.4 







4- 2.9 


4 3.6 


AU-EX 


14.9 





2.8iS: s 8 5 • 10- 2 


4 2.9 


4 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4 3.6 


AL-NI 


11.8 





< 1.8 • 10 3 » 


4 2.9 


4 3.6 


AL-AU-EX-NI 


11.4 







4 2.9 


4 3.6 


AL-AU-NI 


8.1 







4 2.9 


4 3.6 


AU-EX-NI 


7.8 







4 2.9 


4 3.6 


AL- EX-NI 


6.1 







4 2.9 


4 3.6 



continue ^> 



23 



H-10 21 

(Hz^ 1 ) 


config 


(days) 




5j, 5% 

™ b 95%, 


confidence 
yo /o 


interval 

yo/o 






5.0 





3 +2.6 . 1() -3 


4 2.9 


4 3.6 




1 8 





< 1 8 • 10 — 3 


— 29 


— 36 


i nt n 1 


607 8 







— 29 


— 36 


39.81 


AL-NA 


167.1 





1-1 1 1 0.016 
U - 1 -0.017 


4 2.9 


4 3.6 


AL-EX 


100.8 


I) 


- „ -1- n fiO ^ — 9 

1.9Z", y 4 ■ 10 


4 2.9 


4 3.6 


AL-AU 


65.7 





3^3° • 10~ 3 


4 2.9 


4 3.6 


EX-NI 


62.0 





r\ 1 -(-0.019 
u - ± -0.019 


4 2.9 


4 3.6 


EX-NA 


45.5 


1 


„ ^-1-0 04 

U - U -0.041 


4 4.1 


4 5 


AL-EX-NA 


41.4 





_ 


4 2.9 


4 3.6 


AL-AU-EX 


28.3 







4 2.9 


4 3.6 


AL-AU-NA 


19.6 







4 2.9 


4 3.6 


AU-EX 


15.0 





3.2+°;^. 10 


4 2.9 


4 3.6 


AL-AU-EX-NA 


14.2 







4 2.9 


4 3.6 


AL-NI 


11.9 





3j; 2 ' 6 ■ io~ 3 


4 2.9 


4 3.6 


AL-AU-EX-NI 


11.4 







4 2.9 


4 3.6 


AL-AU-NI 


8.2 







4 2.9 


4 3.6 


AU-EX-NI 


7.9 







4 2.9 


4 3.6 




6 2 







— 29 


— 36 




5.1 





o + 2.6 m 1Q -3 


4 2.9 


4 3.6 




1 8 
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